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The fractional quantum Hall effect (FQHE) in the second orbital Landau level at even- 
denominator filling factor 5/2 remains enigmatic and motivates our work. We theoretically consider 
the effect of the quasi-2D nature of the experimental fractional quantum Hall system on a number 
of FQH states (filling factors 1/3, 1/5, and 1/2) in the lowest, second, and third orbital Landau 
levels (LLL, SLL, and TLL, respectively) by calculating the wavefunction overlap, as a function 
of quasi-2D layer thickness, between the exact ground state of a model Hamiltonian and the con- 
sensus variational ansatz wavefunctions, i.e., the Laughlin wavefunction for 1/3 and 1/5 and the 
Moore- Read PfafSan wavefunction for 1/2. Using large numerical overlap as a stability, or FQHE 
robustness, criterion we find that the FQHE does not occur in the TLL (for any quasi-2D layer 
thickness), is the most robust for zero thickness (strict 2D limit) in the LLL for 1/3 and 1/5 and for 
11/5 in the SLL, and is the most robust at finite-thickness (4-5 magnetic lengths) in the SLL for the 
mysterious even-denominator 5/2 state and the presumably more conventional 7/3 state. We do not 
find any FQHE at 1/2 in the LLL for any thickness for the quasi-2D models considered in our work. 
Furthermore, we examine the orbital effects of a non-zero in-plane (parallel) magnetic field finding 
that its application effectively reduces the quasi-2D layer thickness and, therefore, could destroy 
the FQHE at 5/2 and 7/3, while enhancing the FQHE at 11/5, in the SLL. The in-plane field also 
enhances the LLL FQHE states by making the quasi-2D system more purely 2D. The in-plane field 
effects could thus be qualitatively different in the LLL and the SLL by virtue of magneto-orbital 
coupling through the finite thickness effect. Using exact diagonalization on the torus geometry, we 
show the appearance of the characteristic threefold topological degeneracy expected for the Pfaffian 
state. This signature is enhanced by nonzero thickness corroborating our findings from overlap 
calculations. Our results have ramifications for wavefunction engineering, opening the possibility 
of creating an optimal experimental system where the 5/2 FQHE state is more likely described by 
the Moore-Read Pfaffian state with obvious applications in the burgeoning field of fault-tolerant 
topological quantum computing. 

PACS numbers: 73.43.-f, 71.10.Pm 



I. INTRODUCTION 

The fractional quantum Hall effect (FQHE), discovered 
in 1982i, is a quintessentially strongly correlated quan- 
tum phenomenon where electrons in a two-dimensional 
(2D) electron system condense into a strongly interact- 
ing incompressible quantum fluid ground state^i^ with 
fractionally charged quasiparticle excitations which obey 
anyonio^>^, rather than ordinary fermionic or bosonic, 
quantum statistics. The FQHE occurs at low tempera- 
tures in clean (high mobility) 2D semiconductor struc- 
tures under the influence of a strong external magnetic 
field applied normal to the 2D plane of confinement of 
the electron layer. The subject has been studied exten- 
sively during the last 25 years, and reviews can be 
found in the literature. 

In this article we provide a detailed numerical theoret- 
ical study of the orbital Landau level (LL) dependence of 
the FQHE, emphasizing the relative importance of the 
quasi-2D layer width, i.e., the "finite-thickness effect", 
of the electron system transverse to the plane of confine- 
ment in the lowest three LLs. In the non-interacting 2D 
system, taken here to be confined in the x-y plane with 



the magnetic field B along the z-direction which is also 
the direction of confinement with a typical layer width 
of d (finite d corresponds to the realistic quasi-2D sys- 
tem studied in the laboratory and d — corresponds 
to the strictly 2D idealized system often studied theo- 
retically for convenience) , the application of the external 
magnetic field leads to the Landau quantization of elec- 
tronic energy levels given by En = (n -I- l/2)fkuc where 
n = 0, 1, 2, 3, ... is the orbital LL index and lUc — eB/mc 
is the cyclotron frequency defining the equidistant en- 
ergy level separation (i.e., the simple harmonic oscillator 
spectrum) between the 2D LLs. Each LL has a macro- 
scopic degeneracy given by (27rZ^)~^ per unit area where 
I — ^Jhc/eB = yJhJnuJc is the magnetic length (which 
is used as the unit of length throughout). For a given 2D 
electron density of Ns (per unit area) , one has a LL filling 
factor V = iVs/(27r;2)-i = indicating the filling 

of the macroscopically degenerate LLs in the 2D electron 
system, li v < 1, only the lowest (orbital) Landau level 
(LLL), by convention denoted as n = 0, is fractionally 
occupied by electrons. Our discussion, so far, has ne- 
glected electron spin which is equivalent to assuming the 
2D system to be spin polarized (by a sufficiently strong B 
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field, for example). Including spin degeneracy in the pic- 
ture introduces a factor of two since each orbital LL state 
(i.e., n = 0, 1, 2, . . .) can be filled with both up and down 
spins. Incorporating spin in this "trivial" manner (i.e., 
assuming each orbital LL to be occupied sequentially by 
spin up and down electrons), we get that < v < 2, 
2<i^<4, 4<i^<6, and so on, correspond, respec- 
tively, to spin up-down orbital LLs n = (LLL), n = 1 
(the second Landau level, SLL), n = 2 (the third Lan- 
dau level, TLL) . Our goal in this work is to theoretically 
investigate the FQHE in n =0, 1, and 2 and provide 
a critical comparative study of the FQHE in the LLL, 
SLL, and TLL, emphasizing the key role played by the 
quasi-2D layer thickness parameter d (or more precisely 
the dimensionless parameter d/l) in determining the rel- 
ative strength, stability, or importance of various FQH 
states in different orbital LLs. We consider only com- 
pletely spin-polarized (i.e., spinless) FQH states in our 
calculations since the primary fractional states 1/2, 1/3, 
and 1/5 are universally thought to be spin-polarized. 

The motivation of our work stems from the experi- 
mental observation that FQH states are ubiquitous in 
the LLL (about 70 distinct FQH states have been exper- 
imentally observed^'^'^ in the n = LLL with < < 2), 
fairly rare in the SLL (less than 10 FQH states have 
been observedi^iiiii^iHiiiaSiie the n = 1 SLL with 

2 < < 4, and these are much "weaker" than the cor- 
responding LLL FQH states in the sense that the ob- 
servation of the SLL FQHE requires much lower tem- 
peratures and much higher sample mobilities than in 
the corresponding LL situation), and essentially non- 
existent in the TLL (no robust n = 2, with 4 < < 6, 
TLL FQH state has so far been convincingly observed 
exp er iment all y ^ ^ i ^ ^ ) . We establish definitively, in this 
work, that even a qualitative understanding of the higher 
LL FQHE (i.e., for n > 0) must necessarily include the 
finite width effect (in the ^-direction) of the quasi-2D 
electron layer. We point out that the quantitative role of 
the layer width in the FQHE (even in the LLL) has long 
been theoretically know n^^'^^i^°i^^i^^'^'^ . What we show 
in the current work, through the detailed comparison of 
theoretical numerical results in the n = 0,1,2 LLs ob- 
tained on an equal footing within the same model and 
approximation scheme, is that the higher (i.e., n = 1,2) 
LL FQHE has fundamentally different qualitative depen- 
dence on the quasi-2D layer width parameter d/l com- 
pared with the LLL (n = 0) case. In particular, we 
find that a finite value of d/l is essential in establish- 
ing the FQHE in higher LLs whereas in the LLL, finite 
d/l only serves to quantitatively weaken the FQHE! The 
FQH states weaken (strengthen) in n =0 (1) LLs as d/l 
increases from the strictly 2D c? = limit. In the TLL 
(n = 2), we do not find a stable FQHE at all, although 
our limited numerical results show similar trends in n = 1 
and n — 2 LLs. 

The driving stimulus for studying the FQHE physics 
in the n > orbital LLs is, of course, to shed light on the 
enigmatic = 5/2 FQHE, originally observediS exper- 



imentally in 1987 and subsequently confirmed and fur- 
ther studied experimentally repeatedlyiii^i^ over the 
last two decades. The great fundamental significance 
of the 5/2 FQHE cannot be overstated. With the ob- 
vious exceptions of the original discoveries of the (in- 
teger) quantum Hall effect itself and the subsequent 
1/3 FQHEi, the 5/2 FQHE may arguably be the next 
most important experimental discovery in the field. It 
is the only known (so far) exception (for a single 2D 
layer system) to the famous "odd denominator" rule for 
the FQHE, i.e., the FQHE occurs at an odd denomina- 
tor filling factor v = p/q with q an odd integer (and p 
either even or odd) with a concomitant quantization of 
the Hall conductance into a fractionally quantized Hall 
conductivity (Jxy = (p/'z)(e^//i) and a zero (or a deep 
minimum) in the longitudinal conductivity a^x (and also 
in the longitudinal resistivity p^^ — (Jxx\<^xx + c^y]^^)- 
All other (i.e., except for the 5/2 state in the SLL) ob- 
served FQHE states (e.g., v — 1/3,2/5,3/7,... in the 
LLL; V = 7/3, 8/3, 11/5, ... in the SLL) strictly obey the 
odd denominator rule, but the v = b/2 state, with its well 
defined quantized Hall conductance axy = (5/2)(e^/ft.), 
stands in stark contrast to the odd denominator rule. 
We emphasize that the even denominator nature of the 
5/2 FQHE is not only a curious anomaly, it challenges 
our understanding of the FQHE, as developed^ in the 
Laughlin v — 1/q (with q an odd integer) wavcfunction 
and further developed in the Jain composite fermion the- 
ory based wavefunctions^ for the v = p/q type FQH 
states (still with q odd). The odd integer restriction of 
the filling factor denominator q in the "standard" (i.e., 
Laughlin-Jain) FQHE model is inescapable since it arises 
from the Pauli exclusion principle for the electrons. Any 
even-denominator FQHE must thus fall outside the stan- 
dard Laughlin-Jain FQHE paradigm, and must somehow 
correspond to the condensation of bosons (which do not 
obey Pauli principle and therefore allow for even denom- 
inator FQHE) in the v = b/2 incompressible FQH liquid. 

As an aside, it is worthwhile to mention that the 5/2 
state is the only observed even-denominator FQHE, as 
emphasized above, only for single-layer 2D systems. In 
bilayer (or more generally, multilayer) 2D systems, where 
experiments are carried out in two parallel 2D layers sep- 
arated by a barrier (i.e., a double quantum well struc- 
ture), even-denominator (e.g., v ~ 1/2, 1/4) FQH states 
have been observed32i2^iS2ii2£ rather routinely. These bi- 
layer FQH states are theoretically well-understood^ to 
be strong-coupling paired Laughlin states, which were 
postulated by Halperin'^^ some time ago. For example, 
the observed v = 1/2 bilayer FQH state?^:— has been 
shown^ to be the Halperin 331 state'^^, where tightly 
bound pairs of electrons condense into a bosonic Laugh- 
lin state which is allowed to describe even-denominator 
fractions since the Pauli principle does not apply to 
bosons. Such strongly paired Halperin even-denominator 
states are Abelian states, similar to ordinary Laughlin 
FQH states, in contrast to the Moore-Read non- Abelian 
Pfaffian even-denominator stated (see below) which is 
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thought to describe the weakly paired BCS state under- 
lying the 5/2 single-layer even-denominator FQH state. 
We also mention here that no single-layer (as opposed to 
bilayer) z/ = 1/2 FQHE has even been observed exper- 
imentally although we know of no fundamental reason 
ruling out such paired (either strong-coupling or weak- 
coupling) LLL states. 

The leading theoretical candidate for the v — b/2 (and 
its electron- hole counterpart at = 7/2) FQH state is the 
so-called "Pfafhan" wavefunction'^'^ of Moore and Read 
(MR), which has been extensively studied over the last 
17 years^>^^>^^ia^i^d3i445_ fact, there is no 

other proposed viable candidate ground state wavefunc- 
tion for the observed 5/2 FQH state, and as such, a con- 
sensus has emerged that the Moore-Read Pfafhan wave- 
function is the likely description for the enigmatic even- 
denominator 5/2 FQHE. In spite of this consensus, aris- 
ing primarily out of a lack of any other viable alternative, 
there has been a minority viewpoint^i^ questioning the 
vahdity of the MR Pfafhan in describing the 5/2 FQH 
ground state. In addition, a nagging issue of substan- 
tial importance is why there is an incompressible FQH 
state at half filling in the SLL (i.e., = 2 + 1/2 = 5/2), 
but not at = 1/2 in the LLL (or for that matter, at 
half filling in the TLL, i.e. v = A+l/2 ^ %/2) since 
the Pfaffian wavefunction carries no LL index label and 
is presumably an allowed variational description for the 
half-filled LL in any orbital level n = 0, 1, 2, . . .. Impor- 
tant early work^^*^ showed that the Pfaffian wavefunc- 
tion is a rather fragile description of the v = b/2 FQH 
state, and slight variations in the effective interaction be- 
tween the carriers, as, for example, could arise from the 
finite layer width or from changing the orbital LL in- 
dex could, in principle, affect the validity of the Pfaffian 
wavefunction as a suitable description of the v = b/2 
ground state. The "mundane" details of the dependence 
on the effective electron-electron interaction, rather than 
any deep fundamental principle, is, in fact, the reason for 
the SLL 5/2 state to be an incompressible FQHE state 
whereas the corresponding LLL 1/2 state is a compress- 
ible Fermi liquid state. Thus, the orbital LL and the 
finite layer thickness effects (and also perhaps the inter- 
Landau level coupling, an effect we uncritically ignore in 
this work) on the effective electron-electron interaction 
are the key in determining the relative stability of vari- 
ous incompressible FQH states in different orbital LLs. 

In this work, we concentrate on three primary FQH 
states at filling factors 1/2, 1/3, and 1/5 in three or- 
bital LLs n = 0,1,2, and theoretically investigate the 
relative stability of the incompressible quantized Hall 
states at these fillings by calculating, as a function of the 
finite width parameter, the ground state wavefunction 
overlap between the exact (numerically exactly diago- 
nalized) few-particle ground state wavefunction with the 
corresponding candidate variational wavefunction (i.e., 
Laughlin for v =1/3, 1/5 and MR Pfaffian for v ^1/2) 
for the incompressible FQH state. A high (low) overlap 
provides a strong hint that the corresponding realistic 



FQHE state is (is not) described by the corresponding 
variational state (i.e., Laughlin for 1/3, 1/5 and Moore- 
Read for 1/2) in the appropriate orbital LL. Such exact 
diagonalization studies of small systems have been the 
main theoretical tool in learning about the nature of in- 
compressible FQH states ever since the original discovery 
of the FQHE. In particular, the universal acceptance of 
the celebrated Laughlin wavefunction as the appropriate 
description for the observed 1 /3 FQH state in the LLL is 
based almost entirely on the remarkably large (essentially 
unity) overlap between the analytical Laughlin wavefunc- 
tion and the exact small-system numerical wavefunction. 
Similarly, the Jain composite fermion theory based varia- 
tional wavefunctions are thought to be excellent descrip- 
tions for the non-primary (i.e., v = m/{2pm± l)-with 
m > 1 and p an integer-such a,s v = 2/5, 3/7, 4/9, 
etc.) LLL FQH states because of the good overlap be- 
tween the composite fermion wavefunctions and exact nu- 
merical wavefunctions for small systems. Even the MR 
Pfaffian wavefunction is accepted to be the reasonable 
description for the experimental SLL 1/2 state (i.e., at 
1/ = 2 + 1/2 5/2) based simply on the observation, 
made originally by Morf^** and followed up in other subse- 
quent theoretical publications2^i22,, that the Pfaffian has 
good, albeit not spectacular, overlap (~0.8-0.9 for the 
strict 2D system) with the exact small-system numeri- 
cal wavefunction. Conversely, alternative possibilities for 
the 5/2 FQHE state, such as the Halperin 331 state^^ 
or the Haldane-Rezayi spin-singlet hollow-core"'^ state, 
are ruled out^^'^° theoretically essentially entirely on the 
basis of very poor calculated overlap of these candidate 
states with exact small system numerical wavefunctions. 

Another measure of the stability-besides the wavefunc- 
tion overlap-that applies to the MR Pfaffian state, in 
particular, is the appearance of the threefold topological 
degeneracy in the torus geometry that is a signature of 
the Pfaffian state. This degeneracy is a direct signature of 
the non-Abelian nature of the state (and therefore would 
not arise if the state is Abelian such as the Halperin 331 
state). The existence of this degeneracy for a physical 
interaction, i.e., a modified Coulomb interaction, has not 
been earlier reported in the literature. Our finding of the 
appropriate MR degeneracy in the 5/2 numerical state 
precisely where the overlap is maximal is an important 
advance in theoretical understanding. 

We choose the Laughlin (for 1/3 and 1/5) and the Pfaf- 
fian (for 1/2) as our candidate variational ansatz because 
these are the only proposed incompressible FQH states 
at these primary fractional fillings. Our work will, there- 
fore, miss out on discovering any other possible ground 
state wavefunctions (i.e., other than Laughlin for 1/3 and 
1/5 and Moore-Read for 1/2) describing the experimen- 
tal FQH state since we restrict our overlap calculations 
entirely to the Laughlin or Moore-Read candidate states. 
This restriction is, however, not a serious drawback of our 
work since no other candidate wavefunctions exist in the 
literature for the FQH states of our interest. 

Another incentive for our work comes from the fact 
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that existing numerical work^ i^i^i^° , much of it car- 
ried out in the ideahzed 2D hmit, concludes that the 1/3 
FQH state in the SLL (i.e., the i/ = 2 + 1/3 = 7/3 state 
or its electron-hole counterpart v = state) is un- 
likely to be a simple Laughlin state (although in Ref. [1^ 
the deviation from the Laughlin wavefunction stems from 
residual quasiparticle interactions) since the SLL small- 
system exact diagonalization studies give only modest 
wavefunction overlap between the Laughlin wavefunction 
and the exact numerical ground state in the SLL. As 
mentioned above, for the 5/2 FQH state, the calculated 
overlap^^i^^i^ between the MR Pfaffian and the exact 
diagonalization finite system wavefunction is also quite 
modest, much below unity, in the idealized 2D limit. Our 
results would shed light on whether the finite-thickness, 
or width, of the quasi-2D layer enhances the overlap be- 
tween the Laughlin (or Moore-Read) state and the exact 
wavefunction so that the Landau level dependent sta- 
bility of the incompressible FQHE can be discussed in 
terms of a systematic tuning of the Hamiltonian through 
the variation in the electron-electron interaction caused 
by changing the layer width parameter d/l. We find that 
indeed a finite 2D layer width, i.e., a true quasi-2D sys- 
tem, is necessary for stabilizing the FQHE in the SLL, 
at least within the restricted wavefunction space (i.e., 
Laughlin and Moore-Read) we investigate. We note that 
increasing d/l from zero (i.e., the strict 2D ideal limit) is 
equivalent to softening the interaction. 

One rather direct consequence of the quasi-2D finite- 
thickness effect on the FQHE is the orbital coupling of 
an in-plane (i.e., parallel to the 2D layer) appHed mag- 
netic field, in addition to the quantizing perpendicular 
magnetic field. Having both parallel and perpendicular 
magnetic fields is, of course, equivalent to having a tilted 
magnetic field, which is often used in 2D physics to study 
spin polarization effects - the idea being that the paral- 
lel component of the magnetic field couples only to the 
electron spin through the modification of the Zeeman en- 
ergy. This is certainly true in the idealized strictly 2D 
limit where a magnetic field in the 2D plane has no cou- 
pling to the orbital 2D dynamics of the electrons, but 
does couple to the electron spin. In a quasi-2D system, 
however, the situation is more complex since an in-plane 
field could, in principle, couple to the quasi-2D orbital 
dynamics of the electrons through the finite-thickness of 
the electron layer. One effect of the parallel field would, 
for example, be to squeeze the electron layer in the third 
direction, reducing its effective quasi-2D layer width. We 
consider such magneto-orbital coupling effects due to a 
parallel applied field in our work within the context of our 
LL-dependent finite-thickness studies, using a parabolic 
or simple harmonic confinement model. 

A compelling and timely reason for the detailed in- 
vestigation of higher LL FQHE is the recent interest in 
using the v = 5/2 (and possibly v = 12/5, recently ob- 
served by Xia, et al in Ref. [H, and u — 13/5) FQH states 
for fault-tolerant topological quantum computational^ 
using the non-Abelian quasiparticles associated with the 



Pfaffian state. Since the non-Abelian nature of the quasi- 
particle excitations is crucially tied to the specific form 
of the MR Pfaffian wavefunction, it is important that we 
know whether the observed 5 /2 state is really the Moore- 
Read state or not. Recently, serious questions have been 
raised^i^ about whether the Pfaffian is the appropriate 
description for the 5/2 state. Our work, investigating the 
detailed nature of the incompressibility as a function of 
the quasi-2D layer width parameter in higher LLs, thus 
is timely and necessary for further progress in the sub- 
ject of topological quantum computation. There have 
also been several recent experimental investigations of 
the SLL FQHE motivated by topological quantum com- 
putation considerations^^i^i^. 

The rest of the paper is organized as follows: In sec- 
tion \n~K\ we introduce the essential Hamiltonians we 
are studying, namely, the electron-electron interaction 
Hamiltonian and ones that give either the Laughlin or 
Pfaffian states as zero-energy ground states. In sec- 
tion [HB] we describe the various models used to charac- 
terize the quasi-2D nature of the experimental quantum 
confinement, i.e. the "finite-thickness" models for various 
physical systems. The /- and g-functions (defined below) 
are then considered in section III CI as a first attempt to 
understand the effective electron-electron interaction due 
to the quasi-2D nature of the system. Overlaps between 
the exactly diagonalized system for the "realistic" quasi- 
2D system and the Laughlin (fillings 1/3 and 1/5) or 
Pfaffian (filling 1/2) wavefunction as a function of the 
quasi-2D layer width are reported in section UlI Al The 
appearance of the ground state threefold degeneracy sig- 
nature of the MR Pfaffian state on the torus-particularly 
for quasi-2D systems-is investigated in Sec. IIII Bl The 
effect of an in-plane magnetic field and the quasi-2D na- 
ture of the system are considered in section IIII CI In 
section ITVl we discuss connections of our overlap results 
to previous work and, in particular, whether or not the 
physics of the FQHE is adequately captured by the first 
few pseudopotentials. The finite thickness effects on the 
excitation gaps are also briefiy discussed in section |Vl 
Finally, conclusions are given in section IVII with some 
additional discussions. A short letter reporting some of 
our results has recently appeared in the literature^. 



II. MODEL 

A. Hamiltonians: Coulomb, Laughlin, and Pfaffian 

We consider spin polarized electrons entirely confined 
(i.e., no inter-LL coupling) to a Landau level of index 
n = (LLL), 1 (SLL), and 2 (TLL) interacting through a 
pair potential V{rij), where = l^i — '^jl is the distance 
between two electrons (distance is measured in units of 
magnetic length I). Since the electrons are confined to a 
single LL the kinetic energy is a constant, therefore, the 
Hamiltonian for N electrons is taken to be the interaction 
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FIG. 1: (Color online) Vm as a function of m (m odd) for 
the pure Coulomb potential (zero-thickness), i.e., V{k) — 
{e^l/e){l/k), in the LLL (solid circle), the SLL (upward tri- 
angle), and the TLL (downward triangle). Vm are given in 
units of e^/el. The lines are a guide to the eye. 



Hamiltonian 
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H 



(1) 



i<3 



Haldane^i^ showed how this Hamiltonian can be pa- 
rameterized by the relative angular momentum m be- 
tween two electrons, through an expansion in the "pseu- 
dopotential" functions vJi^\ which serve as a complete 
set of basis functions due to angular momentum conser- 
vation: 



N 



N 



E 

=i(odd) 



(2) 



where Pmi'm-ij) is an operator that projects onto the 
states of relative angular momentum rriij — m. Since we 
are considering spin polarized fermions only odd pseu- 
dopotentials are relevant. The Haldanc pscudopotentials 
Vm^"^ for electrons confined to a LL with index n, in the 
planar geometry (as opposed to the spherical geometry) , 
are written as 



dfcfc[L„(fcV2)]'L™(fc2)e-'=V(fc) (3) 



where L„(x) are Laguerre polynomials, and V{k) is the 
Fourier transform of the interaction potential V{r). To 
define our Fourier transform convention we write 



ZTT 



drrJo{kr)V{r) 



(4) 



This parameterization allows all the calculations to be 
done entirely within the Hilbert space of the lowest Lan- 



dau level, i.e., the information about higher LLs is com- 

(n) 

pletely contained within the Vm ■ Note that this simpli- 
fication depends on our neglecting all Landau level mix- 
ing effects, which, along with the assumption of complete 
spin-polarization, is a key and uncritical assumption for 
our theory. 

In a purely 2D system the electron-electron interaction 
is the Coulomb interaction V{r) = (e^/e^)(l/r), where r 
is the distance in the 2D plane between a pair of electrons, 
yielding a Fourier transform of V{k) — (e^Z/e)(l/fc) 
where k is units of l/l. Figure [1] displays the first six 
pscudopotentials as a function of m for the LLL, the 
SLL, and the TLL. Although the differences between 
the pscudopotentials in different LLs are small quanti- 
tatively, the qualitative changes in the system behavior 
can be severe, as discussed below. The main qualitative 
difference in the pscudopotentials Vm^ among the three 

.(1) yi2) 



LLs is that V}"^ > Vl"' > V^"', but V^""' < V. 



.(1) 



(2) 



.(0) 



and V3^°^ < V^'' < V, 

We will denote the celebrated Laughlin^ wavefunction 
at = 1/g to be l^i). It was shown by Haldane^ that 
IvPi,) at filling v = 1/3 is the exact zero-energy ground 
state for a "hard-core" Hamiltonian where 
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H 



(3) 



N 

'E- 



(const.) \^ Pi (rriij ) 



(5) 



This can be obtained from the original Hamiltonian 
(Eq. ID) by setting f/"^ = constant and id"^ = for 
all m > 3; in other words, the Laughlin state at 1/3 
avoids all electron pairs with m < 3 since Vi/Vm — oo 
for m > 3. Further, the Laughlin state at = 1/q (filling 
1/5 for example) is the exact zero-energy ground state 
for a Hamiltonian 



H 



(9) 



(const.) 
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N 



i=i(odd) '<J 



(6) 



Generally, H'j^^ is a Hamiltonian that penalizes two elec- 
trons that have an angular momentum smaller than m. 

The Pfaffian^'^ wavefunction |^p/) is thought to be the 
leading candidate for the FQHE at filling 1/2 in the SLL, 
i.e., v = 5/2, and is known^ to be an exact zero-energy 
ground state for a three-body Hamiltonian Hpf which 
penalizes states where three electrons are in a relative 
angular momentum state smaller than some value. There 
is no known two-body interaction Hamiltonian which has 
the Pfaffian state as the ground eigenstate, so we should 
think of the PfafRan as a variational ansatz for our two- 
body interaction Hamiltonians. 

(We are providing necessarily very brief discussions of 
the Laughlin and Pfaffian Hamiltonians in order to facil- 
itate easier discussions in later sections of this work.) 

In the pure 2D system, where the above discussion 
applies, the electron interaction is just pure Coulomb, 
however, the finite extent of the single-particle electron 
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and (3) the Fang- Howard (FH) variational potential^i^ 
for a heterostructure 



FIG. 2: (Color online) Fourier transform of the finite- 
thickness potentials divided by thickness d, i.e., V{k)/d ver- 
sus kd for the ZDS (solid line), SQ (dashed line), and FH 
(dashed-dotted line) potentials. Also shown is the Coulomb 
potential. 



wavefunction in the direction perpendicular to the plane 
in actual experimental quasi-2D systems will modify the 
electron interaction as discussed below. 



B. Finite-thickness Modeling Potentials 

Depending on the details of the physical systems (i.e., 
quantum wells, heterostructures, etc.), there are many 
possible models for the inclusion of the finite-thickness 
effect with the 2D Coulomb interaction. We do not, how- 
ever, anticipate much qualitative difference among these 
different models since they all lead to the softening of 
the Coulomb potential, with the softening depending cru- 
cially on the largeness of the thickness parameter d/l. In 
particular, the pseudopotentials change in the presence 
of the finite-thickness since, in addition to the magnetic 
length, a new length scale d characterizing the quasi-2D 
thickness becomes relevant, and tor d:s> I, the pseudopo- 
tentials are modified substantially from their d — ideal 
2D Coulomb values. 

We consider the three most extensively used finite- 
thickness potential models, for which we give the Fourier 
transforms, that seek to model the effect of the quasi- 
2D nature of the experimental system, namely, (1) the 
Zhang-Das Sarma (ZDS) potentialSS., which was intro- 
duced specifically to theoretically model finite thickness 
effects on the FQHE, 



VzDsik) 



(7) 



(2) the infinite square-well (SQ) potential^°, which is ap- 
propriate for 2D GaAs quantum well structures, 



(kd)^[{kd)^+4: 



VpHik) = 



eH 9 (24 + 9fcd+ (fcd)2) 
Tsfc (3 + kdf 



(9) 



Potentials (2) and (3) are found by using single-particle 
electron wavefunctions in the 2;-direction of rj{z) = 
^/2ldcos{TTz/d) and 'q{z) = ^27/2d3z exp(-3z/2d), re- 
spectively (these functions are given merely to make our 
definitions of the thickness parameter d, for each model, 
completely clear). In the above, for (1) and (2), d is 
the width of the electron layer (in units of /) in the z- 
direction, and for (3) it parameterizes the layer thickness 
variationally. Obviously, as d — > 0, all of the above finite- 
thickness potentials, describing quasi-2D systems, reduce 
to the pure 2D Coulomb potential V{k) — {e^l/e){l/k). 
Fig. [2] shows the Fourier transforms divided by d as 
functions of kd for each finite-thickness potential used 
clearly indicating the "softening" of the Coulomb poten- 
tial. Note that we keep the background lattice dielectric 
constant e in our definition of the 2D Coulomb interac- 
tion only for the sake of completeness with e^/ (e^) being 
our energy unit (and T the length unit). We note that 
our definition of d as the relevant thickness or width pa- 
rameter for each model of quasi-2D confinement is con- 
tained entirely in Eqs. [7][3]~this is important since later 
we introduce alternative width parameters d' and w. 

As previously mentioned, the purpose of this work is 
not to determine the quantitative accuracy of some par- 
ticular finite-thickness model compared to experimental 
systems. Rather, we are interested in the possible non- 
trivial qualitative changes that can occur when consider- 
ing realistic potentials which are not pure 2D Coulomb. 
We note, however, that the SQ and FH models corre- 
spond to the two most common quasi-2D experimental 
systems (quantum well and heterostructure, respectively) 
whereas the ZDS model, while not corresponding to any 
physical system, is extensively used in FQH theoretical 
studies. 

For the sake of completeness we provide the first six 
pseudopotentials for all three finite-thickness potentials 
as functions of d in the LLL, the SLL, and the TLL shown 
in the left, middle, and right panels of Fig.[3l In all LLs 
shown, the finite-thickness has the effect of reducing (or 
softening) all of the pseudopotentials in a rather "triv- 
ial" way, that is, there is no crossing or non-monotonic 



behavior: V^"' > V^^ > V^'"' > . . . remains for all d. 
It is clear, however, that for the SQ and FH potentials 
the softening as d increases is less severe compared to the 
ZDS potential. (We mention as a cautionary note that al- 
though the same thickness parameter d has been used in 
Figs. [Hand [3] for our three quasi-2D models, the param- 
eter has somewhat different meaning in the three cases 
as can been seen from their strongly different quantita- 
tive effects on the pseudopotential softening in the three 
models as is obvious from Figs.[7][51) Another qualitative 
feature to note is that it is visually difficult to notice any 



.(") 



(8) 
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LLL SLL TLL 

FIG. 3: Pseudopotentials for the LLL (left), the SLL (middle), and the TLL (right). Each panel has three plots each displaying 
the Kn's for the three different finite-thickness modeling potentials, which are labeled: (a) Zhang-Das Sarma (ZDS), (b) infinite 
square- well (SQ), and (c) Fang-Ifoward heterostructure (FIf). The units of the Kn's are as in Fig. [T] 



striking difference between the behavior of the pseudopo- 
tentials, themselves, in different orbital LLs. 

We have actually carried out calculations for a fourth 
model, the parabolic quantum well (PQW) (or the Gaus- 
sian confinement model), which we discuss in Sec. IIII Cl 



C. /- and g-functions 

Only the relative differences in the pseudopotentials 
are important in characterizing the physical nature of 
the FQH ground state. Following Ref. [l^ we form the 
dimensionless /-functions defined through 



(n) _ Yl 



in) 



Vr 



(n) 



(n) 



(10) 



which quantitatively describe how close a given Hamilto- 
nian is to H^^^ that produces the 1/3 Laughlin state as 
the exact ground state. 



From the definition, it is clear that /f 



(») 



(n) 



-1 and 

0. As described in Ref. |l9| the /-function is use- 



ful because for all fm ^=0 for m > 3, and hence, 

any interaction that produces /-functions with this prop- 
erty will be exactly described by the Laughlin state for 
u = 1/3. Further, if the /^-^ are very small for m > 3 
then the exact state will be well approximated by the 
Laughlin state for 1^ = 1 /3, so, by the simple computation 
of the /-functions for some Hamiltonian (or interaction 
potential) one can get an idea of how well the Laughlin 
state will describe the actual ground state. Thus, the /- 
functions far better manifest Laughlin wavefunction-like 



correlations than the pseudopotentials (the Kn"'' func- 
tions) themselves. 

One can generalize this to functions which characterize 
how close a given Hamiltonian is to that produces 
the V — 1/q Laughlin state as the exact ground state. 
To that end we define (^-functions to handle the v — 1/b 
Laughlin state as 



in) 



V, 



in) 



in) 



in) ■ 



(11) 



The g-functions satisfy gj^"'' = — 1 and g^"' — 0. The 

in) 

"hard-core" aspect of a potential is displayed by gm =0 
for all m > 5. One could go further with this procedure 
defining ^.-functions for investigating the "Laughlin-ness" 
of a Hamiltonian at = 1/7, i-functions for — 1/9, etc. 

Figure |4] shows the /-functions as a function of thick- 
ness d for the finite-thickness potentials we have con- 
sidered, namely, the ZDS (top labeled (a)), SQ (middle 
labeled (b)), and FH (bottom labeled (c)) for the n = 
LLL (left panel), the n — 1 SLL (middle panel), and the 
n = 2 TLL (right panel). Note that the /-functions for 
the ZDS potential in the LLL were shown previously in 
Ref. [l^ as well as the /-functions for the FH potential in 
the SLL in Ref. H. For the LLL (left panel of Fig. g]) 

we observe that for zero-thickness d = all /m ^ < 1 and 
it is known (and shown below with other measures) that 
the Laughlin state for 1/3 is a very good approximation 
to the exact state. However, as d is increased all /'s that 
are free to vary {f^^^ through fii^) become larger mono- 
tonically. It is interesting to note that for the SQ and 
FH potentials there is very little increase in the /'s for 
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LLL SLL TLL 

FIG. 4: /-functions for the LLL (left), the SLL (middle), and the TLL (right). Each panel has three plots each displaying 
the /-functions for the three different finite-thickness modeling potentials which are labeled: (a) Zhang-Das Sarma (ZDS), (b) 
infinite square- well (SQ), and (c) Fang- Howard heterostructure (FH). The /-functions, by construction, are dimensionless. 



increasing d compared to the ZDS potential, i.e., ZDS 
overestimates the finite-thickness effect compared with 
the SQ and FH models. This is also apparent in Figs. 
and [3] 

From the calculation of the /m ^ 's, as in Ref.jTol one can 
determine that for some finite-thickness modeling poten- 
tials (e.g., ZDS) the Laughlin state is not a good descrip- 
tion of the physics for d beyond some value whereas for 
other potentials (SQ and FH) the Laughlin state most 
likely remains a good description of the physics for all 
d, but the description becomes progressively poorer with 
increasing d. (This suggestion is further investigated be- 
low by calculating the overlap between the Laughlin state 
and the exact ground state as a function of d.) 

We now investigate /-functions in the SLL: The middle 
panel of Fig. |4] shows /m "* for all potentials considered. 
These functions are very similar to those of the LLL, 
at least qualitatively; they increase essentially monoton- 
ically as a function of d, however, there is a difference. 
Specifically, for d = 0, only /g^^ is below unity while all 
fj^^ through f[]^ are between 1 and 2. This already indi- 
cates that perhaps the Laughlin state will not be a very 
good description of the physics here, and, in fact, this 
has been known for some time^. 

However, there is another property the /-functions re- 
veal as d is increased. The /-functions for the SQ po- 
tential were previously calculated in Ref. 63 and it was 
remarked that in the SLL they do not increase mono- 
tonically the way they do in the LLL. Instead there is a 
weak minimum for intermediate d . That being said, the 
minimum for the SQ potential is weak and /y^'' through 
f[]^ are still greater than unity for all d. Hence, while the 



Laughlin state does become a better description for finite 
d compared to d = it still never becomes essentially 
exact, as in the LLL. In our calculations we observe this 
sort of behavior for the ZDS and FH potentials: an initial 
decrease in the /-functions before a monotonic increase. 

Finally, the right panel of Fig. 2] displays the /- 
functions for the TLL where the behavior follows the 
trend seen from the LLL to the SLL. Namely, /^ 

(2) 

through flj^ are all greater than unity for d = with 

(2) 

fii ~ 3.5. Clearly, the Laughlin state will not be a good 
description of the physics for d = in the TLL. However, 
as d increases the minimum is much more marked and, 
in fact, for SQ and FH potentials at d = /-functions 
are maximum. 

We now consider the ^-functions in order to investigate 
the physics at filling factor 1/5. We would expect the 
Laughlin 1/5 state to accurately represent the physics if 

9m ^ ~ for m > 5. The left and middle panels of Fig. [5] 
display the g-functions for the LLL and SLL, respectively, 
for the three model potentials. It is clear that for d — 
in the LLL and SLL the Laughlin state is likely a good 
representation of the exact state. This is because gin' 
and gm"* are less than 0.5 for m > 5. On the other hand, 
in the TLL (the right panel of Fig. [5]) we observe that, 

f 2) 

for d = 0, 1 > g)n > 0.5, indicating that perhaps the 1/5 
Laughlin state is not a good physical description of the 
exact V = 1/b state in the TLL. 

As a function of the thickness parameter d the g- 
functions behave in much the same way as the /- 
functions, that is, for the ZDS, SQ, and FH potentials 
they increase monotonically as a function of d. The main 
difference between the two is that the local minimum in 
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LLL SLL TLL 

FIG. 5: As in Fig. |3]except for g-functions. 



the g-function is obtained for finite d only in the TLL 
as opposed to the SLL for the /-function (of. Fig. [3]). 
A global minimum is obtained in g for the TLL (right 
panel) for the ZDS potential for finite d while for the SQ 
and FH potentials the minima are obtained for very large 
d. 

It is clear from inspection that one can obtain a rough 
qualitative idea of how well the Laughlin state at 1 /5 will 
represent the physics as a function of d, namely, in the 
LLL the 1/5 Laughlin state will be a good description 
for up to quite large d for the ZDS potential before pre- 
sumably losing out to some other state. The SQ and FH 
potentials, on the other hand, will mostly likely produce 
a state that is very similar to the 1/5 Laughlin for all d. 
A similar prediction is made for the SLL. 

For the TLL, however, the ZDS potential will produce 
a state that is likely not Laughlin-like for small d, a state 
better described by Laughlin for intermediate d, and be- 
comes unlike Laughlin again for large d. The SQ and FH 
potentials are likely to produce a state that is consistently 
unlike Laughlin for all d. 

Our qualitative discussion based on the thickness- 
dependent behavior of /- and gr-functions is consistent 
with experimental findings: The Laughlin FQH state 
is abundant in the LLL, scarce in the SLL, and essen- 
tially non-existent in the TLL. We emphasize that, by 
contrast, no such qualitative discussion is possible with 
respect to the relative abundance of the Moore-Read 
even-denominator state in various LLs since the MR 
state, unlike the Laughlin state, is not an exact ground 
state of any known two-body Hamiltonian, and therefore 
pseudopotential-based functions such as / and g do not 
provide any direct insight into the MR state. 

The /- and g-functions provide a guide to our intu- 
ition and, perhaps, a qualitative understanding. How- 



ever, a way to quantitatively understand the quality of 
the physical description of the Laughlin state is provided 
by calculating overlaps with exact wavefunctions which 
we discuss in the next section. 



III. RESULTS 
A. Overlaps 

A measure of how accurate a variational wavefunction 
(Laughlin or Pfaffian) is compared to the exact wave- 
function is encapsulated in the calculation of the over- 
lap between the two wavefunctions. An overlap of unity 
means the variational state is exact and a vanishing over- 
lap means the variational state is completely unlike the 
exact state, perhaps due to different symmetries (differ- 
ent total angular momentum for example). Although the 
use of overlap calculations in establishing the nature of 
the incompressible FQH states has been a central con- 
ceptual and theoretical tool in FQHE studies, it should 
be emphasized that the overlap calculation has its lim- 
itation since it can only make statements about (nec- 
essarily small) finite systems and specific FQH ansatz 
wavefunctions (e.g., Laughlin, Pfaffian). In spite of these 
limitations the calculation of wavefunction overlap be- 
tween exact numerical wavefunctions for small systems 
with ansatz variational wavefunctions has been a stan- 
dard theoretical FQHE tool for almost 25 years. 

Before we calculate overlaps we make a technical point. 
To calculate properties of the states and diagonalize the 
Hamiltonians we have made use of the spherical ge- 
ometry^2i^iS^we also use the torus geometry later in 
Sec. nil Bl for obtaining the ground state degeneracy at 
5/2 since the ground state degeneracy does not show 
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up in the spherical geometry. The spherical geometry 
is defined by confining N electrons to the surface of a 
sphere with a radial magnetic field produced by a mag- 
netic monopole of strength Q at the sphere center (the 
total flux through the sphere is 2Q{hc/e)). Q is an integer 
or half-integer due to Dirac's quantization condition and 
is related to the radius of the sphere through R = ^/Q. 
The filling factor in a particular LL with LL index n is de- 
fined through its thermodynamic limit v = liniAr^oo ^ I (j 
where g = 2((5 -I- n) -|- 1 is the total LL degeneracy. The 
total angular momentum L is a good quantum number 
and a uniform state is the state with total angular mo- 
mentum L = and is considered incompressible if it has 
a non-zero energy gap between the ground state and the 
low-lying excitation spectra. 

We use the pseudopotentials calculated in the infi- 
nite planar geometry for carrying out our spherical sys- 
tem finite size diagonalization. The reasons for this are 
twofold: (1) it can be argued that the planar pseu- 
dopotentials, since they are the thermodynamic limit of 
the spherical pseudopotentials, better represent the real 
physical 2D system, and (2) it is much more convenient 
when considering finite-thickness modeling potentials to 
use planar pseudopotentials. At any rate, as discussed 
below, it makes very little difference whether we use pla- 
nar or spherical pseudopotentials, and our conclusions 
(although perhaps not the precise values of the overlap 
in each case) are completely independent of this approx- 
imation. We believe that all qualitative conclusions in 
this paper are independent of our planar pseudopoten- 
tial approximation. 



TABLE L Overlap integrals between the exact ground state 
wavefunction using spherical and planar pseudopotentials 
(^*spherel*plane^)' respectively. Also given are the over- 
lap between the Laughlin or Pfaffian wavefunction and the 
exact ground state wavefunction using spherical and planar 
pseudopotentials ({*a|*pia,ne) ^^"^ ^*"l*sphere) *here a de- 
notes either L for Laughlin or Pf for PfafHan as appropriate) . 
This table quantifies the similarities and differences between 
states using planar or spherical pseudopotentials. 



iV 


V 


("^sphere i "^plane ) 


(*-l*plane) 


<*-i*sphere) 
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1/3 


0.998840 


0.992129 


0.996446 




7/3 


0.948005 


0.736947 


0.528481 




13/3 


0.917779 


0.021261 


0.013854 


5 


1/5 


0.999966 


0.996919 


0.997427 




11/5 


0.999974 


0.997886 


0.998198 




21/5 


0.988710 


0.000012 
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1/2 


0.997841 


0.895311 


0.921297 




5/2 


0.968754 


0.963623 


0.867392 




9/2 


0.978191 


0.030311 


0.002384 


10 


1/2 





0.889655 







5/2 


0.972034 


0.934183 


0.837637 




9/2 


0.986467 









For this work we are interested in the FQHE at filling 
factors 1/3, 1/5, an d 1/2 in the LLL, SLL, and TLL. 
As described in Sec. Ill Al all the calculations are done 
in the LLL and all information about higher LLs is en- 
capsulated by the pseudopotentials Vm^ . The relation- 
ship between the LL degeneracy 5 = 2/ + 1 (Z is the 
single-particle angular momentum^^^) and the number of 
electrons N for these states is as follows: = 1/3 has 
2? = 3A^ - 3, = 1/5 has 2l^hN - 5, and v = \l2 has 
21 = 2N — 3. In general, the relation between / and N 
for some LL filling factor ly is 21 = i'~^N + x where x 
is referred to as the "shift" . The shift is non-zero in the 
spherical geometry and a consequence of the finite cur- 
vature of the spherical surface-this is discussed in more 
detail in Sec. UTTB] The {N, I) relation for = 1/2 was 
chosen to be the same as it is for the Pfaffian wavefunc- 
tion^ and the {N,l) for v — 1/3 and 1/5 was chosen 
to be the same as for the Laughlin wavefunction. We 
further note that the FQH states considered in the work 
are not so-called "alias" states, that is, the {N, I) rela- 
tions used only correspond to fillings 1/3, 1/5, or 1/2 
(and 2/3 and 4/5 through particle-hole symmetry). This 
somewhat restricts our choice of particle numbers in our 
finite size diagonalization and will be mentioned again in 
Sec. |V] in relation to the calculation of excitation gaps. 
We believe that the price of restricting N in avoiding any 
aliasing problem is well worth paying in our work because 
our wavefunction overlap calculations become necessarily 
unique since no two distinct FQH states compete at the 
same {N,l) values. 

Confusion about notation can sometimes arise when 
considering FQH states in higher LLs. Here we clarify our 
conventions. In the LLL the 1/3, 1/5, and 1/2 states exist 
at an experimental filling factor of = 1/3, v = 1/5, and 
1/ — 1/2, respectively. In the SLL, the LLL (both spin 
up and down) are filled and inert yielding filling factors 
for 1/3, 1/5, and 1/2 of 1/ = 2 + 1/3 = 7/3, v = 11/5, 
and — 5/2, respectively. Finally, for the TLL, both 
the LLL and SLL (both spin up and down) are filled and 
inert yielding, for 1/3, 1/5, and 1/2, = 4-f 1/3 = 13/3, 
1/ = 21/5, and 9/2. Hence, it should be clear that the 
FQH state at filling 1/5 in the TLL corresponds to an 
experimental = 21/5, for example, and so on. 

Since we are using planar, instead of spherical, pseu- 
dopotentials we quantify the difference between the two. 
Table U provides a number of overlaps for the FQH states 
considered here for the pure 2D {d — 0) Coulomb in- 
teraction. What is clear from these results is that the 
overlap between the exact ground states using either 
the planar (I ^planar)) spherical (I ^sphere)) Pseu- 
dopotentials (Table U third column) are generally high 
(^^spherel^plane) ^ ^■^)- "^^^ major exception is the 
TV = 10 electron system. Here the symmetry of the 
ground state at filling factor 1/2 in the LLL is different, 
on the sphere versus on the plane, yielding a vanishing 
overlap. It is for this reason^ that we do not consider 
the A'^ = 10 electron system in this work. Furthermore, 
the overlap for either the Laughlin or the Pfaffian state 
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FIG. 6: (Color online) Overlaps of either the Laughlin or the Pfaffian wavefunction with the exact ground state wavefunction of 
the pure (d — 0) Coulomb Hamiltonian where Vi (left), V3 (middle), and V5 (right) are varied by 5Vi, SV3, and SV^, respectively. 
Of course, SVi = corresponds to the "Coulomb point" where the interaction is a pure 2D Coulomb interaction. In the above 
figures, the situations considered are u = 1/3 for (A*', I) = (6, 7.5) (dotted line), u — 1/5 for (A'^, I) = (5, 10) (dashed line), and 
V = 1/2 for {N,l) = (8,6.5) (solid line). 



between the exact planar and spherical states is qualita- 
tively similar (Table U columns four and five). One point 
of note is iV = 5 at J/ = 21/5 (1/5 in the TLL); here 
the overlap with the Laughlin state and the exact planar 
state is nearly zero (0.000012) but retains the same sym- 
metry. However, the overlap with the Laughlin state and 
the exact spherical state is exactly zero due to a differ- 
ent symmetry. None of these minor technical issues has 
any bearing on our main goal in this work, which is to 
study the comparative qualitative trends of the quasi-2D 
finite layer thickness effect on the stability of the Laugh- 
lin or the Pfaffian FQH state at = 1/5, 1/3, and 1/2 in 
the lowest three orbital Landau levels. In particular, we 
want to study how the finite thickness of the quasi-2D 
experimental systems affects the comparative stability of 
the Laughlin (for v = 1/3 and 1/5) and the Pfaffian 
— 1/2) state in different {n — 0, 1, 2) orbital Landau 
levels. 

We now return to the calculation of overlaps. In the 
fourth column of Table |T] the calculated overlap between 
the exact state (using planar pseudopotentials) at filling 
factor 1/3, 1/5, and 1/2 and the Laughlin (1/3 and 1/5) 
and the Pfaffian (1/2) states, respectively, for the pure 
2D {d = 0) system in the LLL, the SLL, and the TLL is 
shown for = 5, 6, 8, and 10. In the LLL the overlap at 
1/3 and 1/5 is very high (~ 0.99). For 1/2, however, the 
overlap is not nearly as high (~ 0.9). An overlap of 0.9 
or less is not considered particularly compelling in the 
FQHE and is indicative of, perhaps, different physics. 
In fact, experimentally there is no FQHE observed in 
the LLL at 1/2 to date. In the SLL the story changes. 



For 1/3 the overlap is significantly decreased (~ 0.74) 
while for 1/5 the overlap stays as high as it is in the 
LLL (~ 0.99). The overlap at 1/2 increases in the SLL 
compared to the LLL to a respectable value (~ 0.96). In 
the TLL the overlap at 1/3, 1/5, and 1/2 is essentially 
zero although the symmetry between the two remains 
the same. The most straightforward conclusion follow- 
ing from the d = results shown in Table U is that the 
Laughlin state is the stable FQH state at ly = 1/3 and 
1/5 in the LLL even in the strict 2D limit, but the other 
FQH states may not exist in the ideal 2D limit except for 
1/ = 1/5 in the SLL. 

A theoretical strategy often used in studying the 
FQHE is to vary the first few pseudopotentials away from 
the ideal 2D Coulomb values, in an ad hoc way, to inves- 
tigate whether the overlap between the resulting exact 
states and the Laughlin or Pfaffian states gets better or 
worse. In Fig. [H] we consider a Hamiltonian where we 
have varied Vi, V3, and V5 independently symmetrically 
around the Coulomb point in the LLL (left panel), SLL 
(middle panel) , and TLL (right panel) (the results at the 
Coulomb point are, of course, given in Table [l|. 

Clearly, changes in Vi, V3, V5 that bring the effec- 
tive H closer to ijf ^(ijf ^) increase the overlap between 
the Laughlin wavefunction at 1/3 (1/5) and the exact 
wavefunction, while changes opposite to this decrease 
the overlap. These conditions are obtained, for exam- 
ple, when Vi is increased producing a high overlap for 
1/3 and 1/5 in the LLL and SLL (Fig. [5] left panel (a) 
and (b)). (In the TLL, however, 1/5 has a zero overlap 
due to a symmetry change.) Increasing V3 takes H into 
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^ (5) 

H)^ ' which is evident by the very high overlap between 
the 1/5 LaughUn state and the exact state in LLL, SLL, 
and TLL ((a)-(c) in the middle panel of Fig. Noth- 
ing particularly non-trivial is happening here because we 
know how H is connected to hI^\ i.e., any changes in 

H towards (away from) H^^ makes the Laughlin state a 
better (worse) description for the 1/q FQHE state. 

The overlap variation with changing Vi, V3 , V5 for the 
Pfaflian state {ly = 1/2) apparent in Fig. [S] cannot, how- 
ever, be explained easily since there is no existing two- 
body Hamiltonian, e.g.. Hi for the Laughlin state, for 
which the Pfaflian is an exact eigenstate. It is clear (and 
has been shown previousl y^^'^^ ) that in the LLL and SLL 
changing Vi and V3 in particular ways can produce a state 
with an overlap close to unity for the Pfafhan, i.e.. Fig. [6] 
(a) and (b) in the left and middle panel. It is also inter- 
esting to note that in the SLL the value of 5Vi and SVs 
that gives this high overlap between the Pfaffian and the 
exact state at = 1/2 is very near the Coulomb point. 
Changing V5 (right panel) can also produce a high over- 
lap for the Pfaflian. Qualitatively it behaves similarly 
to the 1/5 Laughlin state as V5 is decreased in the LLL 
and SLL. The Pfaflian never achieves a particularly high 
overlap in the TLL for any values of Vi, V3 , V5 we looked 
at. 

However, the technique of changing Kn's is arbitrar- 
ily artificial and ad hoc, without shedding much light on 
how real experimental quasi-2D systems, where Kn's are 
determined by the layer thickness d, will behave. (In 
fact, all the V„i's change when considering finite thick- 
ness, a fact that is further investigated in Sec. llVI -this 
means that tuning just one specific Vm, e.g., Vi or V3 
or V5, as done in Fig. [6] is purely a theoretical construct 
which is impossible to achieve in real 2D systems.) To 
understand the variation in the states as a function of d 
and their incompressible or compressible nature we care- 
fully define the overlap calculated previously in Ref. [l9l 
which is (\E'i(pj)|\E'j^(£i)). This overlap quantifies exactly 
how similar \^^{d)) is to the Laughlin wavefunction 
(1/ = 1/3 or 1/5) or the Pfaffian wavefunction |5'p/) 
(i^ = 1/2) as a function of the finite layer thickness d. We 
have some intuition (and previous results^^) about how 
the overlap will behave with increasing d after studying 
the /- and 5-functions (cf. Sec. HTCl) . In Figures El [i 
and [5] we report overlaps in the LLL (left panel), SLL 
(middle panel), and the TLL (right panel) of the Laugh- 
lin or Pfaffian wavefunction and the exact wavefunction 
for the finite-thickness models of ZDS, SQ, and FH for 
fillings 1/3, 1/5, and 1/2, respectively. We emphasize 
that the variation in d cannot be described in terms of a 
variation in the value of one (or a few) pseudopotentials 

(K^). 

In Figs. El [51 [21 [ini and [T31 we normalize d defining a 
new width parameter d' such that all the "widths" for 
each potential are defined equivalently. For the SQ and 
FH confinement (as well as the PQW confinement) one 
can calculate w = (z^) — (z)^ which is the variance of 



the wavefunction (i.e., the root-mean-square fluctuation 
in the electron position). Using the values of w for each 
potential, wsq for the SQ and wfh for FH, we normal- 
ize the FH potential to the SQ potential. That is, we 
rescale d in the FH potential to d' — {wfh /wsQ)d — 
(0.57735/0. 180756)d, while d' for the SQ confinement is 
just the original d in that model (since wsq/wsq = 1). 
To rescale the ZDS potential we use a more ad hoc but 
well defined method. Using = 1/2 in the SLL we scale 
the maximum in the overlap between the exact wavefunc- 
tion at d and the Pfaffian wavefunction for the ZDS po- 
tential to be equal to the maximum using the SQ poten- 
tial, i.e., d' — (4.6/1.4)(i. In this way the behavior of the 
overlap as a function of d' for each finite-thickness poten- 
tial is quantitatively similar. Further, in Figs.lTllSll^fTUl 
and[T31we also give the overlaps as functions of w in units 
of magnetic length. With this parameterization one is 
able to distinguish between two regimes of layer thick- 
ness, i.e., w/l < 1 and w/l > 1. Note that our rescaling 
of the width parameter from d to d' (oc w) is a purely non- 
essential book-keeping procedure which makes it explicit 
that, when the quasi-2D width parameter is properly de- 
fined (i.e., d'), then the different quasi-2D models show 
similar quantitative trends in the calculated overlaps as a 
function of layer width. Theoretical descriptions in terms 
of d or d' are completely equivalent-the only advantage 
of using the normalized thickness parameter d' is that 
the calculated overlap is now quantitatively similar in all 
the quasi-2D models we consider. 

For the sake of completeness we provide below the for- 
mulae for w = \/ (z^) — (z)^ and d! for the four models 
(SQ, FH, ZDS, and PQW) respectively in terms of their 
wavefunction parameter, i.e., d, as given in Eqs.EJHKand 
Eq.[T8|) 

Wsq = 0.180756d, 
Wfh — 0.57735(i , 
wpQW = 0.5c? , 



d!. 



SQ 



WSQ 
WSQ 



d — d , 



d'pH = — = 3.194G85d, 



d'. 



SQ 



''ZDS 



WSQ 
WpQW 
WSQ 

4.6 , 
— d = 
1.4 



d = 2.76616(i, 



3.285714d. 



Note that for the ZDS model there is no single particle 
wavefunction in the z-direction, rj{z), that produces the 
effective potential of the form of the ZDS model, hence, 
we cannot define wz ds 1 and d' is calculated as described 
above. 



1. Filling factor 1/3 (Laughlin wavefunction) 

We first concentrate on filling factor 1 /3 (Fig. [7]) • In 
the LLL (left panel) the overlap between the Laughlin 
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FIG. 7: (Color online) Overlap {{^ Ll'^vid))), as a function of thickness d' , between the Laughlin and the exact wavefunctions 
at the fractional filling 1/3 in the n-th LL, i.e., the LLL {ly — 1/3) (left), the SLL {v = 7/3) (middle), and the TLL {v = 13/3) 
(right). The overlap as a function of w is also given on the horizontal axis at the top of each plot. The particular system has 
= 6 electrons at I — 7.5. The different curves represent the different finite-thickness modeling potentials: SQ (solid line), 
FH (dashed line), and ZDS (dashed-dotted line). Note the different scale in each plot of the vertical axis. 



state and the exact state decreases as d' increases. In 
fact, this was previously shown for the ZDS potential 
in Ref. [l^. The SQ and FH potentials also show a de- 
crease in overlap for increasing d' but the change is mild 
compared to the ZDS potential in that the overlap drops 
from near unity only to ~ 0.975-0.96®^. Qualitatively, 
the finite-thickness could eventually destroy the FQHE 
at 1/3 in the LLL, again, as first reported in Ref. ,19.. 
But our results tend to support a scenario where such a 
destruction is continuous as a function of c?', not abrupt 
at any particular value of d' . More specifically, our finite 
system study, as shown in Figs.[7][9l does not indicate the 
occurrence of a thickness-driven sharp transition from an 
incompressible FQH state to a compressible one-any such 
transition seems continuous, although in practice the ex- 
citation gap may become very small for large thicknesses. 

In the SLL (middle panel of Fig. [7|) the overlap starts 
at a modest value 0.74) (cf. Table |J) for d' = and 
an increase in overlap is seen for increasing d' until a 
maximum is obtained for large d' (increasing d' to unre- 
alistic values produces an overlap that approaches zero) . 
We note that the value of d' where the highest overlap is 
obtained corresponds to w// w 1 which, perhaps, more 
clearly shows the effect of the layer thickness, i.e., the 
overlap is seen to decrease as w/l moves away (nega- 
tively or positively) from unity-this qualitative feature is 
again seen when filling factor 1/2 in the SLL is investi- 
gated, cf. Fig. [9l This result could be anticipated from 
the observation that the /-functions for the SLL have a 
local minimum for non-zero d (or scaled d'). However, 
it should be noted that the overlap, while increasing for 



non-zero d' , still only reaches a modest value of approxi- 
mately 0.84-0.92 (depending on the quasi-2D model). 

In the TLL (right panel of Fig. [7]) the overlap never 
reaches a value greater than 0.5 (for any model) indicat- 
ing that the Laughlin state is not a good description of 
the exact state in the TLL. Again, this result is evident 
in the Laughlin-unfriendly nature of the /-functions in 
the TLL. 



2. Filling factor 1/5 (Laughlin wavef unction) 

Next we consider filling factor 1/5 in the LLL (left 
panel) , SLL (middle panel) , and TLL (right panel) shown 
in Fig. [H] From an investigation of the g-functions we 
would expect the 1/5 Laughlin state to have a high over- 
lap with the exact state for a large range of d' in the 
LLL and SLL and that is, in fact, what can been seen 
in Fig. [5] In the LLL and SLL the overlap is greater 
than 0.995 for all models up to large d' (the overlap in 
the LLL and SLL appear nearly identical in their quali- 
tative and quantitative behavior). This result indicates 
that the 1/5 FQH state is as strong in the SLL as it is in 
the LLL, which is consistent with earlier results^'^i^'^. 
Note that nothing interesting happens to the overlap in 
any of the first three LLs a,t w/l k, \. 

The TLL provides a strange scenario for 1/5. Here the 
exact ground state has a different symmetry from the 
Laughlin state at d' = yielding a vanishing overlap. 
However, as d' > the overlap very abruptly becomes 
non-zero and large (~ 0.95). Again, large d' reduces the 
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FIG. 8: (Color online) Same as Fig. [7] except for fractional filling 1/5. The particular system has N = 5 electrons at / — 10. 



overlap severely for the ZDS potential (not shown) while 
only moderately for the SQ and FH potentials (although 
in the reported range of d' , on this scale, the overlap 
appears essentially constant). 

Note that the qualitative behavior of our calculated 
overlap for filling factor 1/3 and 1/5 in the SLL is dif- 
ferent in that for 1/3 increasing d' improves the overlap 
while for 1/5 increasing d' worsens the overlap, as it does 
in the LLL. Based on our overlap calculation we would 
therefore predict a more (less) stable 7/3 (11/5) FQHE 
with increasing quasi-2D layer thickness. 

3. Filling factor 1/2 (Pfaffian wavef unction) 

Lastly, and most importantly, we turn to filling 1/2 
in the LLL (left panel), SLL (middle panel), and TLL 
(right panel) in Fig. [51 considering now the overlap be- 
tween the Moore- Read Pfaffian wavefunction and the ex- 
act numerical wavefunction ai v = 1/2. In the LLL, the 
overlap starts at approximately 0.9 for d' = 0. However, 
unlike the LLL behavior for 1/3 and 1/5 where the over- 
lap monotonically decreases as d' increases, the overlap 
between the Pfaffian and the exact wavefunction here in- 
creases to a weak maximum for some finite d' . The maxi- 
mum overlap for the ZDS, SQ, and FH potentials is mod- 
erate and barely above the d' = value of approximately 
0.9, however. Nevertheless, there is a distinct quasi-2D 
width induced enhancement of the Pfaffian overlap here 
at V = 1/2, not apparent in the corresponding v = 1/3, 
and 1/5 Laughlin states in the LLL where the overlap 
decreases monotonically with increasing layer thickness. 

In the SLL the overlap of the Pfaffian wavefunction 
with the exact wavefunction increases from 0.96 for d' = Q 
to essentially unity for finite d' /I ~ 5. This result sug- 



gests that the finite layer thickness of the quasi-2D exper- 
imental system actually leads to an exact wavefunction 
that is more like the Pfaffian at z/ = 5/2 for finite d' than 
for d' = 0. Continually increasing d' beyond this optimal 
value produces wavefunctions with decreasing overlaps 
(similar to 1/3 and 1/5). This finite width induced stabi- 
lization of the Pfaffian state at v = 5/2 perhaps explains 
the fragility of the observed 5/2 FQHE in experiments. 

In both the LLL and SLL we see that the highest v — 
1/2 Pfaffian overlap occurs for w/l « 1. This is similar 
to what was observed for filling factor 1/3 in the SLL 
(Fig. [T]) but the effect is more pronounced at 1/2. 

Filling factor 1/2 in the TLL behaves much the same 
as for 1/3. The overlap starts below 0.05 for d' — Q 
and achieves a value which is model dependent, between 
0.125-0.25 for some finite d' (for very large d' all overlaps 
eventually approach zero) . From this calculation one 
would not expect to see the FQHE at i/ = 1/2 in the 
TLL. 

The preceding calculations establish that including fi- 
nite layer thickness of the realistic quasi-2D system pro- 
duces non-trivial behavior of the overlap between either 
the Laughlin or Pfaffian wavefunction with the exact 
wavefunction. The physical expectation based on our 
overlap calculation would be that the FQHE in the SLL 
is stronger for the 1/5 and 1/2 state (provided there is 
some finite d') than it is for the 1/3 state. Experimentally 
1/3, 1/5, and 1/2 (and their particle-hole conjugates) are 
all observed in the SLLi^ii^. 

To make a connection between experiment and the- 
ory we consider the recent experimental observations of 
V — 5/2 reported by Choi, et al.— and Pan, et alf^^ 
which were obtained in quantum well structures of width 
30 nms. In Ref. 15 and 16 v — 5/2 was observed at a 
B field of approximately B k, 4.6 T and 5.2 T, respec- 
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FIG. 9: (Color online) Same as Figs. [7] and [8] except for fractional filling 1/2. Hence, the relevant overlaps are between the 
exact ground state wavefunction at some thickness d' {\^v{d,))) and the Pfaffian wavefunction (|^'p/)). The particular system 
has N — 8 electrons at I = 6.5. 




d'Ul 



FIG. 10: Same as Fig. [9] (middle panel, SLL) except only 
the SQ potential is shown. The experimental d' of Refs. 15 
and[lg are shown as a asterisk and solid square, respectively. 
The particular system has N = 8 electrons at I = 6.5. Also 
show is the result for a very recent experiment (Ref. [ssl ) by 
Dean, et al. 



lively. Using the standard formula for magnetic length 
in GaAs/AlGaAs quantmn wells {I ~ 25nm/A/i3[T]) we 
plot in Fig. [10] the value of d' these particular experi- 
ments correspond to on our overlap plot at 1/2 in the 
SLL {ly = 5/2) for the SQ potential. From this compar- 
ison we see that the experimental systems are not opti- 
mized to observe the strongest possible Pfaffian state at 
v = 5/2. Somewhat increasing the value of the quan- 
tum well width so that d'/l (or equivalently d/l) is in- 



creased to the optimal value should give more a stable 
5/2 FQHE! This is further elaborated (and reinforced) 
by studying the ground state topological degeneracy on 
a torus, a defining hallmark of non-Abelian states in the 
next section. 



B. 



Threefold topological degeneracy signature of 
the Pfaffian state on the torus 



Recall (cf. Sec. IIII Al) that in the spherical geometry, 
the signature of an incompressible FQH state is the ex- 
istence of a rotationally symmetric uniform state with 
total angular momentum L = with a finite excitation 
gap to higher energy states. A given FQH state generally 
has a "shift" in the equation relating the number of elec- 
trons N to the total flux 2Q through the finite sphere. 
The MR Pfaffian wavefunction is written on the sphere 



1 



and Vj ~ sin(6'j/2)e*'^^^^ with {9, (/>) being the coordinates 



where the spinor coordinates are Uj = cos{9j/2)e 

being the cc 
The Pfaffian symbol above 



on the surface of a sphere, 
corresponds to 



Pf (Ay) =^e„yl^(i)^(2) ...A,, 



{N-l)cr{N)j 



(13) 



where a are permutations of the N particle indices. It 
is found that the wavefunction in Eq. [T^] requires a flux 
2Q = 2N — 3. While this corresponds to filling factor 
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1/2 in the thermodynamic Umit (A'^ oo) the nontriv- 
ial constant shift of —3 in this relation is characteris- 
tic of the Pfaffian state and arises from the curvature 
of the spherical surface. For example, a possible com- 
peting state at z/ = 1/2 is the composite fermion Fermi 
sea2a with 2Q = 2iV - 2. While the ability to discern 
between competing states via the "shift" may be con- 
sidered an advantage of the spherical geometry, it also 
has the drawback that one cannot directly address the 
competition between different phases without moving to 
a different Hilbert space with a different flux. (In addi- 
tion, complications may also arise from distinct thermo- 
dynamic FQH states occurring at the same shift for given 
{Q,N) values-the so-called 'aliasing' problem.) 

This problem is resolved using the torus geometry 
where the shift is zero and all states are uniquely de- 
fined by the filling factor alone. Thus, exact diagonaliza- 
tion on the torus provides additional information on the 
nature of the ground state. We have performed calcula- 
tions in periodic rectangular domains with unequal sides 
a and b. In the presence of a magnetic field, the standard 
translation operators no longer commute but they do sat- 
isfy the so-called magnetic translation algebra. This non- 
commutation of the standard translation operators pre- 
vents a simple straightforward construction of conserved 
quantum numbers. Haldaneii^ has, however, shown how 
to construct many-body states that have conserved pseu- 
domomenta corresponding to the magnetic translations 
along the two periodicity directions. These pseudomo- 
menta are bi-dimensional {Kx,Ky) and they reside in 
a two-dimensional Brillouin zone containing exactly TVq 
points where TVq is the greatest common divisor of N and 
A^0, where is the number of flux quanta through the 
system (here we denote the total flux as N^jj compared to 
2Q, as in the spherical geometry, to distinguish the ge- 
ometries more readily and because iV^ is more commonly 
used in the torus geometry literature). The pseudomo- 
menta are of the form = ^ixhsja and Ky — i'Kht/b 
with s, i = 0, . . . , A^o- 

On the torus, there is a degeneracy due to the center of 
mass motion given by g at a filling factor p/q. In the con- 
struction of Haldane, the Hamiltonian is block-diagonal 
with exactly q identical blocks-this holds independently 
of the Hamiltonian and, hence, it has nothing to do with 
the physics of the system. We have systematically dis- 
carded this trivial degeneracy in all that follows. The 
construction of conserved quantities has the practical ad- 
vantage that it reduces the size of the Hilbert space in 
which we search-through exact diagonalization-for the 
few low-lying eigenstates. 

At least some of the candidate states for a half filled 
LL have characteristic signatures in the quantum num- 
bers of these low-lying eigenstates. In the case of the 
composite fermion Fermi sea, the effective theory is that 
of quasi- free fermions with an interaction induced mass. 
In a finite system, with discrete energy levels, one expects 
to find closed shell effects which, in addition, should be 
quite sensitive to the aspect ratio of the unit cell. This 



is what is observed"^^ in the LLL at v = 1/2. There are 
also other competing (compressible) phases with broken 
translational symmetry2^i22ii2i which are close in energy- 
the so-called stripe and bubble phases. The stripe phases 
in the LLL are equivalent to charge-density waves with 
unidimensional electronic density modulation. The spec- 
tral signature of such a phase is a set of low-lying states 
with pseudomomenta all related by a single wavevector 
which has the periodicity of the stripe. In the SLL, 
Rezayi and Haldane^^ have numerically shown that such 
a stripe phase is the ground state for electrons interacting 
with the pure Coulomb potential with zero width. 

On the other hand, the Pfaffian state has a very dif- 
ferent spectral signature which can be obtained by trans- 
lating the wavefunction (Eq. I12|) in the torus geometry. 
One of the key ingredients of the Pfaffian is the Laughlin- 
Jastrow correlation factor Wi^j (uiVj — UjVi) which, if the 
coordinates on the torus are z — x + iy, can be written 
as rii<j(-^i ~ -^i)- The correlation factor with the proper 
periodicity is given by 9i{zi — Zj|t) with r = ib/a {b/a 
is the aspect ratio of the rectangle) and 9i the first Ja- 
cobi theta function. This allows a construction^''' of the 
standard Laughlin state for filling factor 1/m (with an 
additional factor to treat the center of mass motion). 

This same construction does not work for the Pfaffian 
state. Indeed, the presence of a denominator in the for- 
mula (Eq. [T2)) invalidates the periodicity properties if one 
simply introduces 9i{zi — Zj\T) factors. The correct pro- 
cedure^^'^^ involves the four Jacobi theta functions 0i, 
62, 03, and 64 through the following substitution. 



This construction gives three ground states instead of 
only one in the spherical geometry. This remarkable de- 
generacy is of topological origin. 

If one considers a physical Hamiltonian whose ground 
state is adiabatically connected to the Pfaffian state 
one expects to find the threefold degenerate multi- 
plet whose quantum numbers on the torus may be de- 
duced from Eq. 1141 In a finite system, there will be 
a splitting of the multiplet of states, and it is only 
in the thermodynamic limit that these states will be- 
come truly degenerate. For electrons at filling 1/2 the 
three Pfaffian ground states have pseudomomenta K — 
(0, A^o/2), (A^o/2, 0), iNo/2, No/2) (in units of 2T:h/a and 
2'KTi/h) and this set of quantum numbers clearly differ- 
entiates the Pfaffian from the other ground states men- 
tioned above. Hence, to probe for this quasi-degeneracy 
one has to use a rectangular unit cell since, in that 
case only, the degeneracy is non-trivial-square or hexago- 
nal cells have additional discrete (geometric) symmetries 
leading to the equivalence of some or all Brillouin zone 
points potentially masking the non-Abclian topological 
degeneracy. 

To fully understand the properties of the Pfaffian state, 
one further has to take into account the fact that it 
breaks particle-hole (p-h) symmetr y ''^1^^ , either due to 
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FIG. 11: Exact energy in units of e^/el (only the low en- 
ergy sector is being shown) versus pseudomomentum K= 
^/{a/b)Ki + {b/a)K^ for iV = 16 electrons, using the torus 
geometry, interacting via the SLL Coulomb Hamiltonian (zero 
width). The pseudomomenta Kx and Ky are given in units 
of 2irh/a and 2nh/b, respectively. Panel (a) corresponds to 
an aspect ratio of the rectangular unit cell equal to 0.99 while 
(b) has an aspect ratio of 0.75. 
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explicit p-h symmetry breaking terms present in real sys- 
tems (e.g.. Landau level mixing or coupling) or due to a 
spontaneous p-h symmetry breaking in theories neglect- 
ing LL mixing^*'. If we consider the wavefunction written 
on the sphere, its flux 2Q = 2N — 3 leads automatically 
to the p-h conjugate state (the so-called anti-Pfaffian) at 
2Q = 2N + 1, which, since it exists at a different flux for 
the same N, precludes any mixing between the states. 
On the torus these two states will mix, for finite size 
systems, which should lead to a lower energy symmet- 
ric Pfafhan- Anti-Pfaffian combination and a higher en- 
ergy antisymmetric combination. This reasoning applies 
to each of the three states with differing K's expected 
on the torus, so we expect to find three doublets if we 
are, in fact, dealing with a system well described by the 
Pfaffian state, corresponding to a non-trivial topological 
degeneracy of 6 in the thermodynamic limit. 

To investigate the spectral signature of the Pfaffian 
we have performed exact diagonalizations on the torus 
from TV = 10 to 16 electrons using the SLL Coulomb 
interaction with finite thickness modeled by the three 
quasi-2D confinement models used before-FH, SQ and 
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FIG. 12: Same as Fig. [11] except here we have included fi- 
nite thickness effects using the SQ (a), FH (b), and ZDS (c) 
potentials, respectively. The aspect ratio is 0.75. 



ZDS. 

For small systems, there is no obvious threefold de- 
generacy at zero width while the threefold degeneracy is 
clearly seen when finite thickness is included^. For the 
largest system we have been able to study i.e. = 16, 
the picture becomes clearer. For A'^ = 16 electrons inter- 
acting via the pure Coulomb potential in the SLL we find 
that there is a threefold quasidegenerate set of ground 
states with quantum numbers of the Pfaffian for aspect 
ratio 0.99: see Fig. [TlTa). However, we do not observe 
candidate, higher-lying, states forming their p-h doublets 
that are expected from the higher energy combination of 
Pfaffian- Anti-Pfaffian states. If we tune the aspect ratio 
from 0.99 to 0.75, then the three quasidegenerate ground 
states are still obtained but there is a less clear separation 
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FIG. 13: Same as Fig. [12] except the aspect ratio has been 
tuned to 0.99. 



from the higher-lying states, see Fig. [TlTbV 

If we now consider systems with a finite width taken 
to be the value of maximum overlap (d/l ~ 4-5), as de- 
termined from the calculations on the sphere, then we 
find the Pfaffian signature (threefold degeneracy and p- 
h partner states) is qualitatively enhanced. At an as- 
pect ratio equal to 0.75 the three finite thickness mod- 
els considered all lead to the threefold quasidegeneracy 
for the three states with the correct quantum numbers: 
see Figs. IT^ aVfc). Now it is plausible to identify the 
p-h partners for each of the states with the same quan- 
tum numbers predicted by the Moore-Read wavefunc- 
tion slightly higher in energy than the ones forming the 
threefold quasidegenerate ground state manifold. If the 
aspect ratio, with finite thickness included, is changed 
from 0.75 to 0.99 we find that things change somewhat. 



see Figs. rTST aVfc): while there is still a threefold quaside- 
generacy of the ground state, the doublet structure is no 
longer clear. The reason for this sensitivity to the as- 
pect ratio is not known and deserves further study. We 
emphasize, however, that even if the enhancement of the 
Pfaffian signature is restricted to some range of aspect 
ratio, its presence is clearly enhanced by finite thickness. 

The results presented in this subsection, taken with 
those in Sec. IllI Al provide a satisfying picture where 
finite thickness effects produce a ground state that is 
described by the MR Pfaffian wavefunction. The fact 
that two different geometries produce the same qualita- 
tive conclusion using two different, complementary, sig- 
natures is quite convincing. On the sphere (Sec. IIII Al 
and Figs.[7][5|) we find that the overlap between the Pfaf- 
fian and the exact finite-system numerical wavefunction 
for the ly — 5/2 state is enhanced substantially as the 
quasi-2D thickness parameter is increased. On the torus 
(Sec. HITB] and Figs. [TW^ we find that the expected 
non-Abelian topological degeneracy, a characteristic sig- 
nature for both Pfaffian/anti-Pfaffian states, shows up 
precisely where the wavefunction overlap is large. 



C. In-Plane magnetic field effects: Overlaps 

The consequence of the application of an in-plane mag- 
netic field to the ^ = 5/2 FQH state is a very important 
question. Experiments have shown^L^ that the FQHE 
at z/ = 5/2 is suppressed with increasing in-plane compo- 
nent of the magnetic field. Originally, this was thought 
to point towards a spin-unpolarized FQH state at 5/2, 
since traditionally the in-plane field is assumed to cou- 
ple only to the spin degrees of freedom, and increas- 
ing the in-plane field is supposed to enhance the spin- 
polarization of the system. However, subsequent theo- 
retical workSi seemed to settle the debate regarding the 
spin-polarization of the 5/2 state, and the 5/2 state is 
now considered to be spin-polarized. A question, there- 
fore, arises about the suppression of the 5/2 FQHE in- 
duced by the finite in-plane field since spin-polarization 
presumably cannot play a role in a spin-polarized FQH 
state. It is speculated that the suppression arises from 
the in-plane field induced orbital effects leading to mod- 
ifications of the pseudopotentials V,n ■ 

The in-plane field serves to squeeze the single-particle 
electron wavefunction in the direction perpendicular to 
the plane. Thus, the application of an in-plane field 
serves to effectively decrease the width of the quantum 
well. From the previous section (i.e., section Ull Ap it is 
clear that such a reduction in the effective width could 
lead to the destruction of the FQHE at v = 5/2 and 
7/3 whereas it could seemingly strengthen the FQHE at 
ly = 11/5. 

We now explore this in detail. Here we find that a 
parabolic confinement model is convenient, i.e., the z- 
direction confinement potential is V{z) = (2?i^/TO(io)z^ 
where do is the thickness of the quantum well. (In this 
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FIG. 14: (Color online) Overlap between the exact ground 
state wavefunction using a parabolic confinement at filling 
factor 1/3 (dashed line) and 1/2 (solid line) and the Laughlin 
and Pfafiian wavefunction, respectively, as a function of d' , in 
the SLL (top panel) and the LLL (bottom panel). 



section we assume that the do of the actual system is 
such that the overlap between the exact ground state 
and the Pfaffian wavefunction is maximum, that is, do is 
a constant and not a variable, cf. section UlI A 31 ) The 
single-particle ground state electron wavefunction in the 
z-direction is thus a Gaussian 



(15) 



We then apply an in-plane magnetic field of strength Bi 
with the vector potential Ai = (BiZ, 0, 0) and now make 
a simplifying assumption, namely, we ignore the cross 
terms (i.e., we consider the limit Bi <^ B where B is the 
perpendicular magnetic field strength) thus arriving at a 
slight modification, from the original quantum confine- 
ment, of the Schrodingcr equation for ri{z) as 



2^ + ^(^) + ^4F 



viz) = Eri{z) 



(16) 



where If = fic/eBi. This equation has a solution of the 
same form as Eq. 1151 i.e., 



ry(z) 



(17) 



with l/ct^ = l/dl + l/Alf. This imphes that the quasi-2D 
layer thickness d in the presence of i?i 7^ is less than the 



original thickness do > d, i.e., the 2D layer is squeezed 
by B,. 

The potential used when calculating the planar pseu- 
dopotentials is 



V{q) 



e2 erfc(gd/2)e(9'*)'/4 



(18) 



which includes the effects of a parabolic quantum well 
confinement and an in-plane magnetic field through 
the definition of d above. Note that the effective Si- 
dependent apparent width d decreases as Bi increases, 
squeezing the layer width parameter in the presence of 
the in-plane field. 

The parabolic confinement is qualitatively similar to 
the confinement models used previously (SQ, FH, ZDS) 
and, in fact, quantitatively almost identical to the Fang- 
Howard confinement (see Ref. ,37.1 . In Fig. [14] we show 
the calculated overlap between the exact ground state 
wavefunction for 1/3 and 1/2 in the LLL and SLL with 
the Laughlin and Pfafhan wavefunction, respectively, as 
a function of d', to convince the reader that the difference 
between using a parabolic confinement and the Fang- 
Howard confinement is very small provided d' is defined 
as it was in section UlI AI (here d' = (0.5/0. 180756)d scal- 
ing it again with respect to the square- well confinement). 

Obviously, V{k)/d vs. fed, Vm vs. d, and /- and g- 
functions vs. d (cf. Figs. [1] [31 21 and [5]) are very similar 
to the results shown for the Fang-Howard confinement 
and are, therefore, not shown. 

To investigate the effect of a finite in-plane magnetic 
field we find it convenient to plot overlap as a function 
of the dimensionless variable A — {Bi/B) between the 
range < A < 2. This variable is r elated to the effective 
thickness d as A = 2^1/d'' — 1/dg, where do is chosen to 
be the value of thickness (without any in-plane magnetic 
field) at which the overlap is the highest, i.e., dg/l ~ 5. 

Finally, in Fig. [TH we plot the overlap between the 
exact ground state wavefunction for a parabolic confine- 
ment at a finite in-plane magnetic field Bi and the Pfaf- 
fian wavefunction for 1/2 in the SLL (i^ = 5/2), the 
Laughlin wavefunction at 1/3 in the SLL {v — 7/3), and 
the Laughlin wavefunction at 1/3 in the LLL {v — 1/3), 
respectively. The application of an in-plane magnetic 
field (i.e., Bi,X ^ 0) causes the exact state at = 5/2 or 
7/3 to become more unlike the Pfaffian or Laughlin state, 
respectively, since the overlap goes down monotonically 
with increasing Bi/B, and hence, the in-plane field could 
eventually destroy the FQH at v = 5/2 and 7/3 by ef- 
fectively enhancing the confinement, making the system 
more two-dimensional. On the other hand, the in-plane 
field causes the Laughlin state to become a better candi- 
date for the exact state at 1 /3 since it makes the system 
more 2D (this is also true for i/ = 1/5 in the LLL and 
SLL, which is not shown). Note that all of our in-plane 
field results follow directly from our findings in Sec. lIII A] 
where we showed that increasing the layer width stabi- 
lizes the Pfaffian 5/2 and the Laughlin 7/3 states whereas 
it destabilizes the Laughlin 1/3 state. In this section, we 




FIG. 15: Overlap between Pfaflian and the exact ground state wavelunctions at 5/2 (left), the Laughlin and the exact wave- 
functions at 7/3 (middle), and the Laughlin and the exact wavelunctions at 1/3 (right) as a function of A. 



explicitly show that an applied parallel field, by squeez- 
ing the quasi-2D layer width, could suppress the 5/2 and 
7/3 states but strengthen the 1/3, 1/5, and 11/5 states. 

As important word of caution is, however, in order 
with respect to our in-plane field results presented in 
this section. We have only considered one particular 
aspect of the applied in-plane field, namely, the quasi- 
2D confinement effect through Eqs. [TSlfTSl There are 
several additional effects induced by the in-plane field, 
which, although not considered in our work, may very 
well be important in practice. Two obvious spin ef- 
fects of the in-plane field, neglected in our work be- 
cause we are only considering completely spin-polarized 
states, are the Zeeman coupling induced spin polariza- 
tion of the ground (and possibly excited) states. More 
subtle magneto-orbital effects neglected in our work are 
the in-plane field induced orbital anisotropy in the 2D 
plane and the subband-Landau level coupling induced 
enhanced scattering due to the finite in-plane fields, and 
the possible in-plane field induced stabilization of com- 
peting compressible states (e.g., stripe or bubble phases) 
which may have lower ground state energies than the 
incompressible FQHE states under our exclusive consid- 
eration in this work. A full consideration of all possible 
effects of the in-plane magnetic field is well beyond the 
scope of this paper, where we have concentrated entirely 
on the wavefunction overlap effect for FQH states arising 
from the wavefunction squeezing (i.e., d < do) by the ap- 
plied in-plane field. In real experimental situations, some 
of these neglected effects may very well be significant or 
perhaps even dominant. 

IV. CONNECTION TO PRIOR WORK 

The main question that we address in this section is 
the following. Do the first few Haldane pseudopotentials 
determine the FQHE physics in the SLL? This has been 
earlier studied in the literature, and therefore this sec- 
tion connects our work to the existing work of Moi&^ 
and Rezayi and Haldane'^^, which have been influential 



in theoretical studies of the 5/2 state during the last 10 
years. 

As mentioned in the Introduction (Sec.|T]), early impor- 
tant work was done by Mor^ and Rezayi and Haldane^ 
regarding the Pfaffian description of the 5/2 FQHE. In 
particular, the overlap between the exact ground state of 
the SLL Coulomb potential, where the first pscudopoten- 

tial Vi^"^ was varied around the Coulomb point, and the 
Pfaffian was calculated as a function of this variation of 
the pseudopotential, i.e., SVi^\ Rezayi-Haldane utilized 
the torus geometry varying both V-^^"^ and V^^\ and as 
such, their results are not directly comparable to the re- 
sults presented here. Morf, however, varied V-l^^^ using 
the spherical geometry and his results are directly appli- 
cable to ours. In fact. Fig. 3(b) in Ref. is equivalent 
to our Fig. ^h) (middle plot)-of course, we have also 
included overlaps with respective Laughlin states as well 
as considered two (LLL and TLL) other LLs. A general 
result of Refs. [H and [35l was that a positive Sv}^^ en- 
hanced the overlap. However, as discussed in Sec. IIII Al 
for the Pfaffian, the change in overlap via the variation 
in the short range pseudopotentials is not easily moti- 
vated since the Pfaffian is not an exact eigenstate of a 
Vi-only two-body Hamiltonian, and, in fact, an increase 
in v}^'' (leaving all other Vm^'s constant) is physically 
untenable, i.e., there is no experimental or physical way 
one can effect an increase only in Vj^^^ in real systems. 

Here we show that the first few pseudopotentials do 
not by themselves determine the physics. In Figs. Uni fTTl 
and [TSl we show the calculated overlap between the Pfaf- 
fian wavefunction and the exact ground state wavefunc- 
tion of the SLL Coulomb Hamiltonian where {v}^\ V^^''), 
{Vi^\v^^^), and {V^^\v^''''') are allowed to vary away 
from their original SLL values by SVm^ for N — 8 elec- 
trons and 2{Q + n) = 13, i.e., for the 5/2 state. Each 
plot has left and right panels. The left panel is a color 
contour plot of the overlap as a function of the change 
in pseudopotentials. The right panel displays the same 
contour plot but with only contour lines at values of the 
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FIG. 16: (Color online) Contour plot of the overlap between 
the PfafBan wavefunction and the 5/2 ground state wavefunc- 
tion of a SLL Coulomb Hamiltonian, where V\'^' and V^^^ have 
been varied from their pure Coulombic values, as a function 
of the variations ^Vj^^' and 5V,}^\ The system considered is 
the N = 8 electron system shown in Fig. [9l The left panel 
is the color contour plot while the right panel contains only 
contour lines at values of the overlap equal to 0.5, 0.8, 0.9, 
0.95, and 0.99. Also shown in the left panel is the path traced 
out in the {V^^\ V^^'^) plane of the finite thickness of the SQ 
potential from d/l = to 10. 
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FIG. 17: (Color online) Same as Fig. 1161 expect here we con- 
sider variations in v}^^ and V^^'' . 

overlap equal to 0.5, 0.8, 0.9, 0.95, and 0.99. Also on this 
panel is a series of square points that show the overlap 
for the SQ potential at different values of d/l (we use d 
here since d' — d for the SQ by definition) from zero to 
d/l = 10 in unit steps. Finite thickness changes all the 
V^'^'s but we are only showing the path traced out in 
the V^^'^-V^'^'^ or V^^'^-V^^'^ or V^-V^^'^ pseudopotential 
space of the contour plot. 

Two main observations can be made in the parameters 



FIG. 18: (Color online) Same as Fig. 1161 expect here we con- 
sider variations in V^^^ and V^^\ 

space of {VI'\VP), (f/^\ F^^^), and {V^^\^ V^). First, 
there are always regimes where the overlap is above 0.99. 
That this occurs in the {yi^\v^^'^) space is not particu- 
larly surprising considering the previous work of Rezayi- 
Haldane'^^ (although this is shown here for the first time 
in the spherical geometry). However, the other two con- 
tour plots in the {V^^\v^^^) and {V^^\v^^'') spaces go 
against the conventional wisdom. In particular. Fig. 1181 
shows that a large overlap with the PfafRan can be ob- 
tained by leaving V^*-^^ constant and varying V^^^ and 

V^'^^ only. Hence, it is clear that the first two pseudopo- 
tentials do not dominate the physics of the 5/2 FQHE. 

The second observation that can be gleaned is that, 
since the values of all the pseudopotentials change upon 
including finite thickness effects, one cannot parameter- 
ize finite thickness corrections in terms of only two pseu- 
dopotentials, be they the first and third, first and fifth, 
third and fifth, or any other combinations of two. This 
effect is clearly shown in the right panel of each contour 
plot where the path of the overlap, as a function of finite 
width, in the {V^'\V^% {Vi'\,Vp), and {y^^\v^'^) 
plane, respectively, does not coincide with a high over- 
lap region shown in the contour plot when only chang- 
ing two pseudopotentials. Note that the squares of the 
finite thickness overlap results are shaded according to 
the color contour plot on the left panel (as well as being 
more clearly visible in Fig. [9]). 

We have thus demonstrated that the investigation of 
the overlap as a function of changes (or ratios) between 
two pseudopotentials does not elucidate the physics of 
the 5/2 state and produces ambiguous results. In partic- 
ular, in the inset of Fig. 5 of Rezayi-Haldane2^ the path 
of the finite thickness is plotted in the {Vi'^\ V^^'^) plane 
and shown to cross a boundary between a striped state 
and a paired FQHE state. However, from the above con- 
siderations, one cannot conclude that the effects of fi- 
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nite thickness drive the system across boundaries in this 
way. Specifically, finite thickness corrections involve vari- 
ations in all the pseudopotentials from the ideal Coulomb 
point, and tuning one or even a few pseudopotentials does 
not, under any circumstances, mimic the finite thickness 
effect. Tuning one or two pseudopotentials around the 
Coulomb point in order to study the stability of the 5/2 
FQHE state with respect to the Pfafhan is, therefore, 
somewhat misleading in our opinion, particularly since, 
unlike the LLL Laughlin states, there is no theoretical 
two-body Hamiltonian for which the Pfaffian is an exact 
ground state. 
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The FQHE transport activation gap can be readily 
measured experimentally and is an extremely relevant 
quantity that characterizes the incompressibility of a 
FQH state. In this section we discuss the effects of fi- 
nite thickness on the excitation gaps for the FQHE at 
v = 7/3 and 5/2. Usually the experimental activation 
gap and the theoretical excitation gap are considered to 
be the same although this may not be necessarily true 
in the presence of disorder. It is, however, well-accepted 
that larger FQHE gaps imply stronger FQH states asso- 
ciated with larger FQHE excitation energies. 

For finite size systems there are a few ways to calcu- 
late the gap and we only consider the gap in the spheri- 
cal geometry. The gap is considered to be the energy of 
a well separated quasiparticle-quasihole pair where the 
initial ground state is incompressible. Hence, one only 
considers the excitation gap if the ground state of the N 
particle system at flux Q is a uniform state with total 
angular momentum L = 0. If this is the case then the 
excitation gap can be calculated as the energy of a quasi- 
particle E^''^'^^^^ and a quasihole E^^^^^^^^ where E^g'^^ 
is the ground state energy of a system of N electrons at 
flux Q-hence, the quasihole state has Q + 1/2 while the 
quasiparticle has Q — 1/2. With these energies the gap 

is given as A E^^'^^^^^^ + E^'^'^^^^ - 2E\^\ Note that 
this definition involves only calculations of ground state 
energies albeit at different flux values. 

Another way to find the gap is to calculate the energy 
of the N electron system at flux Q as a function of angular 
momentum L. The ground state (again only if the state 
is uniform and thus incompressible) will have L — Q, 
and a branch of low energy excited states at different 
i 7^ 0. The gap is then given as the energy in the long 
wavelength limit which corresponds to the lowest energy 
excitation with L = N for the state at ly = 7/3 and 
L = N/2 for the, presumably, paired state at = 5/2. 
In other words, 

A = (iex) - E'-'^^ = 0) , (19) 

where L^x = N {N/2) for fiUing factor 7/3 (5/2). This 
definition involves obtaining the low-energy spectra of 



FIG. 19: Energy gap A for f — 5/2 using the SQ potential as 
a function of d' /I (left) and w/l (right). The solid line is the 
gap in units of e^/el while the dashed line is the gap in units 

o^e^/e^/p + d''^ (left) and e'^/e^P + (right). The system 
considered has N — 8 electrons. 



the system at a given Q. This is the method we use in 
this work to investigate the gap. The reason is that, as 
discussed by Morf^^ and in Ref. ^5^, the first method, de- 
scribed above, leads to some ambiguity. This is because 
for v = 5/2 and = 8 electrons the quasiparticle state 
is aliased to a FQHE state with filling factor 2/3 in the 
SLL. Thus, what is being assumed as the quasiparticle 
energy from the 5/2 FQH state may in fact be a 2/3 
filled incompressible FQH state. Since we are also not 
calculating the gap in the thermodynamic limit, it is not 
particularly important which method we use since we are 
only interested in broad qualitative features. See Refs.lsvl 
and 44 for more thorough numerical investigations of the 
energy gaps in the FQHE. We emphasize that our re- 
sults are only qualitative and should not be compared 
quantitatively with experimental activation gaps. 

In Figs, [in] and we plot the calculated gap for 
I' = 5/2 and 7/3 for the SQ model for systems with 
N — 8 and 6 electrons, respectively. The left panel shows 
the gap versus the width d' and the right panel shows the 
gap versus w, both in units of I. The solid fine is the gap 
A in units of e^/el, and we see that for both filling factors 
the gap (in units of /el) decreases with increasing width 
(in units of ^). This is not surprising since the scale of 
the energy itself is being reduced as the Coulomb energy 
is suppressed below /el in the presence of finite thick- 
ness, e.g., see Fig. [3] where the pseudopotentials them- 
selves are shown to decrease with increasing width. In 
fact, this introduces a subtle point when considering the 
theoretical energy gap including finite thickness effects. 
For every value of d! we are, essentially, considering a 
different Hamiltonian, so it is not quite appropriate to 
think of the gap as a function of the well width. This is 
most clear when one considers that the energy scale itself 
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FIG. 20: Same as Fig. 1191 except for p — 7/3 and a system of 
= 6 electrons. 



is changing with d' . To readily incorporate the effect of 
a varying energy scale with varying d' we should rescale 
the gap energy. The dashed lines in the left (right) panels 

of Figs, [inland HOI show the gap in units of e'^/e\/P + d'^ 
{e^/e^/P + w"^). In these units the gap is seen to increase 
with increasing width, i.e., in reduced energy units, the 
excitation gap is enhanced with increasing width. 

We note that two alternative, and not necessarily 
equivalent, signatures exist for the theoretical charac- 
terization of the stability or the strength or the robust- 
ness of FQH incompressible states. These are (1) the 
overlap between the ground state finite-size numerical 
wavefunction and a candidate incompressible state, and, 
(2) the FQH excitation gap calculated directly numeri- 
cally. The relationship between these two signatures of 
incompressibility is not obvious at all, and both defini- 
tions have their problems. The definition of the excita- 
tion gap is ambiguous (particularly due to the aliasing 
problem) since the two ways of defining it, as discussed 
above (using the ground state energy at different Q and 
the excitation spectra at the same Q), are inequivalent. 
Also, a finite size system always has discrete energy lev- 
els, and thus always has a gap. In addition, the exci- 
tation energy A, being an energy, is by definition not 
dimensionless, and therefore is somewhat problematic as 
a signature for the ground state compressibility. (For 
reasons discussed above, we believe that the widespread 
practice of expressing A in units of e^/eZ, the so-called 
Coulomb energy, is arbitrary, and (A/(e^/eZ)) may not 
necessarily be a faithful representation of the stability of 
ground state incompressibility, particularly in situations, 
e.g., 5/2 state, where finite quasi-2D width is essential 
for the existence of the FQHE.) On the other hand, the 
signature of the overlap with a candidate analytical in- 
compressible state (i.e., a variational ansatz such as the 
Laughlin or the MR wavefunction) has several intrinsic 



problems: (1) such an analysis is necessarily limited by 
the constraint of the variational ansatz, i.e., if a differ- 
ent (unknown) incompressible state describes the ground 
state better, the overlap calculation would miss that; (2) 
finite-size effects inherent in small system calculations, 
where a large overlap may turn out to be a misleading 
finite-size artifact; (3) the competition with nearby com- 
pressible states is not captured. 

It is our contention that when one investigates the 
FQHE theoretically at some filling factor v it is most 
illuminating to first determine the nature of the state re- 
sponsible for the FQHE. This is usually done via an over- 
lap calculation with some trial wavefunction (Laughlin, 
composite fermion, MR Pfaffian, etc.) such as was done 
in Sec. IIII Al Once a satisfying identification is made for 
the FQH ground state using the wavefunction overlap 
signature then experimentally relevant quantities, such 
as the excitation gap, can be calculated and compared 
with experiment. In fact, this is essentially the historical 
record of events following the discovery of the = 1/3 
FQHE and subsequent explanation by Laughlin. This is 
also the strategy we have followed in the current work. 



VI. CONCLUSION 

In this work, we investigate the relative stability of 
the FQHE at the most important primary filling factors 
J/ =1/2, 1/3, and 1/5 in the three lowest orbital Landau 
levels n = (LLL), 1 (SLL), and 2 (TLL) by calculat- 
ing, as a function of the quasi-2D layer width parameter, 
the wavefunction overlap between the directly diagonal- 
ized exact many-body wavefunction (for small spherical 
systems) with the corresponding candidate theoretical 
ansatz for incompressible states, namely, the Laughlin 
wavefunction for v = \/?> and 1/5 and the MR Pfaffian 
wavefunction for v — 1/2. Rather surprisingly, we find 
the layer thickness parameter, often neglected in theoret- 
ical studies of the FQHE mostly carried out in the ideal- 
ized zero-thickness strict 2D limit, to be a key parameter 
in stabilizing the incompressible states in the SLL. In par- 
ticular, we find that the SLL states at = 2 -I- 1/3 = 7/3 
and 2 + 1/2 = 5/2 have larger exact overlap with the 
corresponding theoretical incompressible states for finite 
values of the layer thickness parameter d'; typically, the 
overlap is maximum (of the order of unity, in fact) for 
d' /I ~4-5 in the SLL whereas in the LLL the overlap 
decreases monotonically with increasing d' /I, i.e., the in- 
compressibility is the strongest at d' = (note that the 
width parameter d' cx c? is the "normalized" width pa- 
rameter, and d' — d only for the square-well confinement 
potential). Furthermore, we also find that finite width 
is essential to the observation of the threefold degener- 
ate ground state signature of the MR Pfaffian state in 
the torus geometry. In particular, the expected three- 
fold non-Abelian topological degeneracy for j/ = 5/2 only 
shows up on the torus for d' ~ 4-5Z, where the corre- 
sponding overlap with the MR wavefunction is optimal. 
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In the TLL, the overlap is always (i.e., for any d! /l) very 
small, indicating the essentially generic absence of incom- 
pressible FQH states in n > 1 orbital LLs. 

Our finding of the "absence" (i.e., small overlap) of the 
Laughlin v — l/Z state in the d! — Q limit in the SLL 
is consistent with earlier theoretical work carried out in 
the ideal d' = limiiM'^'^^. The experimental FQHE 
observationai^iHiii at = 7/3 and 8/3 occur, of course, 
at finite values of d'/Z, and the experimental values of d! /I 
in the real samples is d' /I ^^2-5, which is in agreement 
with our theoretical finding. Our work indicates that 
tuning the quasi-2D width parameter significantly (either 
far above or below the optimal value where the overlap is 
maximum) should destroy incompressibility at v = 7/3 
(or 8/3). 

The most important motivation for studying the SLL 
FQHE is to understand the nature of the enigmatic 5/2 
FQH state, the only even-denominator FQH state ever 
observed in single-layer 2D systems. What does our work 
imply for the stability of the anomalous FQHE at half- 
filling? We briefly discuss below the qualitative answer 
to this question. 

Our work clearly shows the importance of the finite 
quasi-2D layer width in producing the incompressible 
FQH state at v = 5/2, assuming the state to be de- 
scribed by the MR Pfaffian wavefunction which is the 
only available spin-polarized candidate wavefunction for 
the 5/2 state. Consistent with earlier theoretical work in 
the literature^ii^i^S, we find the overlap between the ex- 
act (small-system) many-body wavefunction at i/ = 5/2 
to have a modest overlap (~0.8-0.9) with the Pfaffian 
wavefunction, which increases monotonically to a large 
overlap of almost unity as d' /I increases. For larger d' /I, 
the overlap decreases again. Thus the behavior of the 
5/2 state as a function of finite width is similar to the 
2 -\- 1/3 = 7/3 Laughlin state-both are most stable at a 
finite layer width in contrast to the LLL Laughlin states 
at 1/3 and 1/5, which are most stable (i.e., maximum 
overlap ^ 1) in the strict 2D {d' /I — 0) limit. 

The behavior of the LLL 1/2 state is interesting in 
this context. The overlap with the Pfaffian here shows a 
very slight increase as a function oid' /I before decreasing 
again similarly to the SLL 5 /2 state except that the Pfaf- 
fian overlap never approaches unity for the LLL 1/2 state, 
indicating that an incompressible = 1/2 FQH state, at 
least one that is well-described by the MR Pfaffian, is un- 
likely to occur in the LLL. We emphasize that there is no 
compelling fundamental reason for the LLL 1 /2 state not 
manifesting a = 1/2 FQHE, it is merely absent in the 
reasonable parameter space of the Coulomb interaction 
for realistic experimental systems. It is conceivable that 
a clever tuning of the interaction Hamiltonian far from 
the strict 2D Coulomb point will stabilize a LLL v = 1/2 
FQHE corresponding to the MR Pfaffian state. 

The v = \/b state in the SLL (i.e., z/ = 2-^1/5 = 11/5) 
is very similar to the LLL FQH states in its dependence 
on the layer width. The overlap with the Laughlin state 
decreases monotonically with increasing width parame- 



ter, implying that both the LLL and SLL 1/5 FQHE 
would be the strongest in the strict 2D ideal limit for zero 
layer thickness, similar to the situation for the = 1/3 
LLL state. 

Our work indicates that, in principle, wavefunction en- 
gineering should be possible to enhance the FQHE at 
V — and 7/3 in the SLL by increasing the quasi- 
2D width parameter d' /I. A trivial way of enhancing 
d' /I is, of course, to decrease I (at fixed d') by increas- 
ing the magnetic field. Since increasing the magnetic 
field B increases the interaction energy ~ /I a/B, 
it is obvious that this would strengthen the FQHE triv- 
ially (the energy gap scales as e^// thus increasing the 
gap). But increasing magnetic field while keeping v fixed 
requires a proportional increase in the carrier density 
which is problematic (and may lead to the occupancy of 
the second subband, reducing the sample mobility sub- 
stantially). Therefore, we suggest that d' should also be 
increased to produce enhanced stabilization of the SLL 
FQHE. The fact that different models of quasi-2D con- 
finement, e.g., heterostructure (Fang-Howard), square- 
well, parabolic quantum confinement, and the Zhang- 
Das Sarma model, all give very similar qualitative and 
quantitative variation of our calculated overlap with the 
width parameter indicates that our qualitative conclusion 
is model independent, and applies to all quasi-2D physi- 
cal systems where the FQHE is experimentally studied. 

As a direct consequence of the layer width dependence 
of the FQHE found in this work, we also consider the ex- 
perimentally important situation*^'^^ of the effect of an 
in-plane magnetic field, B^, applied parallel to the quasi- 
2D layer on the FQHE. Assuming the system to be com- 
pletely spin-polarized, as we do throughout this work, 
and neglecting all disorder effects (i.e., ignoring for the 
moment that Bi could cause additional scatteringS^i^i by 
inducing coupling between in-plane and transverse dy- 
namics), the only effect of Bi is to modify the quasi-2D 
layer width d' by shrinking it to d'^{Bi) < d'^{Bi — 0). 
This in-plane magnetic field induced modification of the 
quasi-2D layer width arises from the magneto-orbital cou- 
pling between in-plane and transverse electron dynamics 
in a finite-width system, which is strictly absent in the 
ideal 2D limit where, in the absence of spin (i.e., for a 
completely spin-polarized system), there can be no cou- 
pling between the in-plane and the transverse motion. 
We show that in a quasi-2D system, the in-plane field 
Bi could destabilize the SLL 5/2 Pfaffian state by de- 
creasing the effective layer width through the reduction 
of the overlap between the exact wavefunction and the 
Pfaffian wavefunction. This follows naturally from our 
finding that in the SLL, the reduction of the quasi-2D 
layer width always reduces the overlap of the 5/2 state 
and since Bi reduces the effective value of the layer width 
it would naturally suppress the overlap at ly = 5/2. The 
same is also true for v = 7/3, but not for v = 11/5 in 
the SLL or for =1/3 or 1/5 in the LLL since the lat- 
ter three FQH states are the most robust (i.e., maximum 
overlap) at the smallest value of d' . We therefore predict 
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that, to the extent the 1/2, 1/3, and 1/5 FQH states are 
completely spin-polarized, the application of a parallel 
magnetic field is likely to weaken the v =5/2 and 7/3 
FQH states and strengthen the v =11/5, 1/3, and 1/5 
FQH states. Of course, the physical effects neglected in 
our approximation (e.g., spin. Landau level coupling, dis- 
order) could play important roles in actual experiments, 
but we suggest careful experiments in an applied in-plane 
field to validate (or falsify) our prediction of the weak- 
ening (strengthening) of the 5/2, 7/3 (11/5) FQHE in 
the SLL. We mention that the quantitative effect of the 
in-plane field depends strongly on li/l and U/dl where 
li = [chj eBi)^!'^ and / and d! are the 2D magnetic length 
and the Si = value of the width, respectively. Since 
lijl and lijd! both should be not too large for the orbital 
influence of the in-plane field to be observable, it may 
not be easy to observe the predicted effect in the LLL 
(or even for the 11/5 SLL) since the typical "/" is rather 
small for these cases. 

Before concluding we summarize the large number of 
approximations we have made in our theory: (1) We have 
assumed a spin-polarized system throughout, and there- 
fore if spin is playing an important dynamical role in 
any of the fractional states we study, then our results 
would not be particularly useful in understanding the cor- 
responding experimental observations. (2) We have ne- 
glected the Landau level coupling throughout our work, 
and it is possible that the LL coupling plays a role in 
the SLL FQHE**^ (3) We have neglected all disorder 
effects-in particular, the application of the in-plane par- 
allel field may introduce additional "disorder" by open- 
ing new channels of scattering (e.g., intersubband scat- 
tering) through the coupling of the in-plane and trans- 
verse dynamics. This will then serve to further weaken 
the 5/2 and 7/3 FQHE in the presence of the in-plane 
field. We expect such scattering effects to be very weak 
in the li ^ d! and li ^ I regime of experimental interest. 
(4) We have uncritically assumed that the "strength" of 
the incompressibility at a given filling factor (or equiva- 
lently, the robustness of the FQHE at a given fraction) 
is determined by the calculated overlap between the ex- 
actly diagonalized numerical wavefunction for small sys- 
tem sizes and the corresponding candidate theoretically 
postulated FQH state, i.e., the Laughlin wavefunction 
for V — 1/5, 1/3 (in all LLs) and the MR Pfaffian wave- 
function for V — 1/2 (in all LLs). (5) Related to the last 
point, our work will completely fail if the actual state de- 
scribing the FQHE at a particular fraction is qualitatively 
different from our assumed candidate wavefunctions. For 
example, one cannot, as a matter of principle, rule out 
the possibility, unlikely as it may seem, that the observed 
V — hl2 FQH state is the MR Pfaffian state for finite 
d! jl ~4-5, as we find, which continuously and adiabat- 
ically goes over to some other unknown incompressible 
state at lower values of d'/L Such a continuous crossover 
implies that the observed 5/2 FQHE will remain strong 
and robust for all values of d', but our calculated overlap 
with just the Pfaffian decreases as d! jl decreases. We 



believe such a scenario to be extremely unlikely, partic- 
ularly since there are no other proposed wavefunctions 
for 5/2, but we cannot rule this out on purely theoretical 
grounds. (The same also applies for the 7/3 state where, 
if anything, such a scenario of two distinct states, Laugh- 
lin for layer d! jl and "something else" for smaller d! jl, 
seems even more unlikely.) (6) We have used approxi- 
mate models for considering the finite layer width in the 
quasi-2D system. However, the fact that four distinct 
models of quasi-2D systems produce essentially identical 
results and conclusions indicate that our results should 
have qualitative and semi-quantitative validity in real 
samples. Also, earlier theoretical worb^ii^^i^ indicates 
that more sophisticated models of quasi-2D confinement 
do not lead to appreciable changes in the FQHE numer- 
ical results compared with the relatively simple models 
used in our work. (7) We have utilized the standard 
spherical geometry for our finite size diagonalization us- 
ing rather modest system sizes (number of electrons be- 
tween 5 and 10 depending on the filling factor). We be- 
lieve that the small system size of our exact diagonaliza- 
tion study is not a problem since we are not attempting 
any quantitative estimates of the excitation gap (or other 
experimental quantities), but are interested in the qual- 
itative dependence of the overlap as a function of the 
layer width in different LLs. One reason for our use of 
relatively modest system size is, of course, computational 
ease (since we are producing a very large amount of nu- 
merical data: four different models of confinement, and 
three different LLs, for three different FQH states, i.e., 
36 sets of diagonalization done as a function of the width 
parameter d!\). Yet a second important reason for our 
choice of system sizes is to avoid the well-known "alias- 
ing" problem on the sphere where two distinct fractions 
occur together exactly at the same parameter values. We 
have only chosen system sizes where the aliasing issue 
does not arise. The fact that we study three different 
orbital LLs on an equal footing to compare the relative 
qualitative stability of the 1/3, 1/5 and 1/2 states gives 
us confidence in our numerically obtained trends (as a 
function of d!) and conclusions, if not the precise num- 
bers. 

Finally, we emphasize an important qualitative find- 
ing of our work. We find that, in contrast to the Laugh- 
lin states (e.g., v = 1/3) in the LLL where increasing 
the first pseudopotential Vi compared with the Coulomb 
value necessarily strengthens the FQHE in a theoretically 
and physically meaningful manner, no such simple pseu- 
dopotential adjustments (either a single pseudopotential, 
e.g., Vi, V3, V5 or arbitrary combinations of them, e.g., 
{Vi,V3), (^1,1^5), (VsiV's)) make sense for studying the 
FQHE stability in higher (SLL or TLL) orbital Landau 
levels. This is because of the theoretical fact that the 1/q 
Laughlin state becomes and exact LLL eigenstate of the 
effective two-body Hamiltonian h'^^ (see Eq. |6]) where 
all pseudopotentials m > q are taken to be zero, e.g., 
the Laughlin 1/3 state is just an exact LLL eigenstate of 
the Coulomb interaction if Vi — > cx). This mathematical 
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simplicity provides an adiabatic connection between the 
Laughlin state and the exact i' — 1/3 FQH state even 
in the reaUstic system as long as Vi is not too "small". 
This mathematical simplicity is, however, completely lost 
in higher orbital LLs where the exact state cannot sim- 
ply be written down in this manner. In particular, the 
V = 5/2 even-denominator state in the SLL, if it is in- 
deed the Moore-Read Pfafhan state, is not an eigenstate 
of any (even a completely unrealistic) two-body Hamilto- 
nian. As such, trying to understand the nature of the 5/2 
state by varying the few lower pseudopotentials (e.g., Vi, 
V3, etc.) around the Coulomb point is not mathemat- 
ically (or physically) well-motivated since no two-body 
Hamiltonian exists with the MR state as its exact ground 
state. We have therefore adopted the physically moti- 
vated approach in this work by working directly with 
the realistic Hamiltonian including the quasi-2D finite 
thickness effects (instead of using unrealistic arbitrary 
variations in the pseudopotentials). Varying the finite 
quasi-2D layer thickness involves changes in all the pseu- 
dopotentials which cannot be mimicked by varying any 
one or two pseudopotentials only. We find that this phys- 



ical approach leads to a rather unexpected finding: While 
in the LLL, finite quasi-2D thickness always weakens the 
FQHE, in the SLL, the finite thickness may actually sta- 
bilize the FQHE, for example, bX v — 5/2 and 7/3 filling 
factor. 

In summary, we have theoretically considered the effect 
of orbital dynamics on the stability of the primary FQH 
states at 1/2, 1/3, and 1/5 filling factors by calculating 
the wavefunction overlap (and the topological degeneracy 
expected for the MR Pfaffian state) as a function of the 
quasi-2D layer width {d!) in n = (LLL), 1 (SLL), and 
2 (TLL) orbital LLs, finding that the FQHE does not 
occur in the TLL (for any value of the quasi-2D width 
parameter), is the most robust at d' = in the LLL and 
for the 11/5 state in the SLL, and is the most robust at 
d' /I --4-5 in the SLL for the 5/2 and 7/3 state. We also 
do not find any v ~ 1/2 FQHE in the LLL for any value 
of the layer width. 
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